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Introduction 


1. The 2-body problem is the name usually given to the 
problem of the motion of a system of many particles 
attracting each other according to Newton’s law of gravita- 
tion. This is the classical problem of mathematical nat- 
ural science, the significance of which goes far beyond the 
limits of its astronomical applications. 

The 2-body problem has been the main topic of celestial 
mechanics from the time of its inception as a science. 
Now that many of the problems of celestial mechanics 
have become part of geophysics, the central position of 
the n-body problem has been further strengthened. 

The fundamental dynamical problem for a system of 2 
gravitating bodies is the investigation and predetermina- 
tion of the changes in position and velocity that the 
particles undergo as the time varies. However, this 
is a very complex non-linear problem whose solution has 
not been possible under the present-day status of math- 
ematical analysis. 

It is natural that first to be developed were the prac- 
tical aspects of the »-body problem, arising in connection 
with the study of the motion of the planets, comets, and 
satellites. These problems were solved by the method 
of successive approximations. However, the mathematical 
difficulties entailed by these problems were so great, that 
the convergence of the approximations was not proved 
nor an estimate of the error obtained. Although the 
computational methods were not rigorously justified, their 
applicability, nevertheless, has been verified by practical 
usage which has shown that there exists an amazing agree- 
ment between theory and observation in most of the cases 
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considered. Despite their mathematical imperfection, the 
practical methods of celestial mechanics yield an accurate 
picture of the motion of the celestial bodies over a limited 
interval of time. 

2. Today, the hope for achieving some success in the 
study of the 2-body problem is centered mainly on the 
use of qualitative methods of investigation. 

A qualitative investigation of the 2-body problem 
usually consists in studying an incomplete set of parameters 
which, though not defining the state of the system uniquely, 
nevertheless, characterizes it in certain respects. Clearly, 
the problem of the predetermination of the values of an 
incomplete set of parameters cannot be well-posed, and 
we must confine ourselves to ascertaining the general laws 
by which these parameters change. This restricted for- 
mulation of the problem, characteristic of qualitative 
investigations, has the advantage that a knowledge of all 
the integrals of the motion is not needed for its solution, 
and significant results can be obtained by using only the 
known integrals, or by analyzing directly the equations 
of motion and their various analytical consequences. 

However, in choosing an incomplete set of parameters, 
it is important that we choose a representative set capable 
of reflecting the essential characteristics of the motion, 
and which leads to meaningful ideas about the dynamical 
system. To the set of such parameters belong, first of all, 
the mutual distances between the bodies of the system. 
It is therefore not accidental, that the search for the most 
general properties of the distances between the bodies as 
the time tends to infinity (positively or negatively) has 
occupied a prominent place in the qualitative investiga- 
tions of celestial mechanics. This problem is often called 
the “problem of final motions’’. 
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The first general and very simple results on the final 
motion of 2 gravitating bodies were obtained by Jacobi 
in 1842. However, after Jacobi’s work, attempts were 
made to study only the special case of three bodies. A 
systematic investigation of the final motion of three 
gravitating bodies was begun about thirty years ago by 
J. Chazy. For this case, he attempted to classify the 
typesof motions possible and toinvestigate their connection. 

The results of Chazy proved to be very interesting, 
although one of his most important ones, the impossibility 
of capture, is incorrect. After O.Yu. Schmidt showed, 
with the help of the techniques of numerical integration, 
that capture is possible, the need arose to carry on a more 
detailed study of the final motion in the three-body 
problem. It was necessary to obtain effective qualitative 
theorems that give sufficient conditions for the occurrence 
of specific types of motions in the form of restrictions on 
the initial data. A general approach to the solution of 
this problem was presented by the author in the years 
1948-1951, and he obtained the very first criteria for the 
occurrence of hyperbolic and hyperbolic-elliptic type mo- 
tions in the three-body problem. 

Research in this direction was continued by G.A. 
Merman, K.A. Sitnikov, V.F. Proskurin, N.G. Kochina, 
G. Ye. Khrapovitskaya, and O.A. Sizova. Thus, by 
today, a considerable amount of scientific work on this 
question has been accumulated (see the bibliography at 
the end of the book). 

We stated above that after Jacobi’s investigations, no 
attempts were made to obtain new results on the character 
of the final motion in the general case of an arbitrary 
number of bodies. The mathematical difficulties of this 
problem have apparently caused researchers to stay away 
from it. Nevertheless, we believe that progress in the 
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general case is possible, although it seems it may not be 
rapid. 

Below, we are going to consider some qualitative 
methods of analyzing the 2-body problem, as well as some 
of the results obtained by means of these methods. Of 
course, the totality of these results does not, by far, give 
a complete qualitative picture of the motion of n grav- 
itating bodies, and its clarification will require still further 
work. 


3. Let us give a brief outline of the content of the book. 


In the first chapter, we give the equations and general 
integrals of the n-body problem, and we study the simplest 
theorems on the final motion due to Jacobi. 

In the second chapter, we consider means of applying 
the method of dimensional analysis to the n-body problem. 
As far as we know, dimensional analysis has not been used 
before in the investigation of the final motion; though it 
does not yield definitive results in this area, it is very use- 
ful in carrying out a preliminary analysis of the problem. 

In the third chapter, we present our “method of con- 
tinuous induction,” and we consider some applications of 
it in which the final motion in the 2-body problem is 
analyzed. This method allows us to obtain effective 
qualitative results, namely, it allows us to formulate suff- 
cient conditions for the occurrence of certain types of 
final motions in the form of conditions on the initial state 
of the dynamical system. 


The fourth chapter is devoted to the method of 
invariant measure. The application of this method to 
the many-body problem has necessitated working out a 
number of theorems on the measure theory of dynamical 
systems. These are presented at the beginning of the 
chapter. At the end of the chapter, we prove some very 
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general theorems on the motion of a system of gravitating 
bodies using the method of invariant measure. 


In the fifth chapter, an attempt is made to analyze 
some cases of the evolution of a system of m gravitating 
bodies on the basis of celestial mechanics. Here, we shall 
be concerned with the processes which are of cosmogonical 
interest, and which are accompanied by the conversion of 
mechanical energy into non-mechanical forms. 
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CHAPTER 1 


Equations and General Integrals 
of the n-Body Problem. 
Simplest Theorems Concerning the Final Motion. 


1. We shall consider a system of n particles P;, Pe, 
...,£, with masses 721, Me,..., M, attracting each other 
according to Newton’s law of gravitation. We shall study 
the relative motion of these particles. At first, we consider 
the motion of P,, Po, ..., P, with respect to the center 
of mass of the system. 


Let us take a set of rectangular coordinates with origin 
fixed at the center of mass of the system and with coordinate 
axes having fixed directions. Let x;, y;,2; be the coordi- 
nates of the point P;,7; the distance of P; from the origin, 
and 7;; the distance between P; and P;. We shall use a 
system of units for which the universal gravitational 
constant is 1. 


Let us set 
Vics Qing 
oa (1.1) 


iév 


where the summation is taken over all pairs of distinct 
bodies; then the equations of motion are given by the 
following 3” differential equations of second order: 
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dt UU, ys __ OU dtz,__ OU 4.2 
at — On Mm, -7- aa —~ Oys ? m, ‘aa at? Fi ( ) 


¢=21,8, 2. 0 
(In this chapter, we shall only consider the equations of 


motion and their integrals.) These equations may also be 
written as 67 equations of first order 


dx, dy; __, dx, 

Ms Me 

U ‘ : 1.3 
dz, 9U dy; au dz; au u2) 


|| a | 
oad dx,’ * dt dy (dt ox 
é=1, 2, va oe B 


A first integral of these equations is the familiar law 
of conservation of energy 


LS m,(z3 + y" +2, ‘y= U -+-H. (1.4) 


=f 


where # is the constant of integration. 
Besides this, there are the six linear momentum in- 
tegrals. In our choice of coordinates they are given by 


3 m.z. =), ym, =0, Vmz,=0, 


(==) = 


ym =0, Sih V m,2, 50. 


(=) i= 


(1.5) 


There are still the three area (moment of momentum) 
integrals which, in terms of the coordinates under con- 
sideration, are 
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i 
p2 M,(2Y,—- YX) =%, 
= 


»2 m; (y 2, = Z.Y,) — Coy (1.6) 


‘=1 


% 
2 m, (2,2, a 22.) = 'Css 


Here ¢;, ¢2, and ¢; are constants of integration. 


These ten integrals exhaust all of the known essentially 
independent integrals of the given dynamical system. 


2. Another theoretically useful way of describing the 
relative motion of 1 bodies can be given by introducing 
Jacobi coordinates. We first change the notation of the 
preceding section into one that is better suited to the 
application of Jacobi coordinates. 


Let the particles in question now be denoted by Pa, 
P,,..., P,-1 and their respective masses by mo, m, ..., 
My)» 


Let &1, 71, £1, be the coordinates of P; with respect to 
axes with origin at the point Po, &, 2, 2 the coordinates 
of P, with respect to axes with origin at the center of mass 
of the particles Py and Py etc., and finally, let &,—1, 
Nw~t, (n—1 be the coordinates of the point P,—1 with 
respect to axes with origin at the center of mass of the 
particles Po, Pi,..., P,~2. The coordinate axes in all 
of the systems are assumed to be parallel. Let p; denote 
the distance of P; from the origin of the coordinate system 
in which the motion of this particle is being described. As 
before, we let 7;; denote the distance between P; and P;, 
and we again use a system of units in which the universal 
gravitational constant is 1. 
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The equations of motion in Jacobi coordinates are then 
given by 
di, = aU d2rg aU dl, __ ou 
a TT -T ~ “OR, ? Ms 3 di2 On ° + Bs diz ~~ al, d (1.7) 


61, 2, .2 0) wml 
where #41, /2,..., My—1 are the reduced masses defined by 


mom, 


py =—=.-s ——— ee 
: mMg+m, ’ 


(mg + my) mo 
My + Mm, + Me * 


hy = 


(my + m+... + Tent) Mn—1 
oar a Mo + my +... + May 


Equations (1.7) may be expressed as 6(m—1) equations of 
the first order: 


dF, __ be dr; Pree avs , 
dt =, a ae ae =, 
dt, au dy, aU de; au (1.8) 
Mt de — 0b PR ae Oy: Pi aE OG 
¢=1, 2, ..., n—1 


In Jacobi coordinates, the energy integral becomes 


an) 


2 SHG tC v+H, — as) 
Gan 


where #7 is the constant of integration. 
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The momentum integrals are satisfied identically in 
Jacobi coordinates, and the moment of momentum in- 
tegrals take on the form 


n=] 


Py HF, — 18) =%,, 


n—! 
; ; 1.10 
2 B,(no,— C7) = es, ( 


u—] 


Py OF, — oo) =e, 


where ¢1, €2, and ¢s are constants of integration. 


3. Finally, we give one more way of describing the 
relative motion of m gravitating bodies Po, Pi,..., Px—1 
using a coordinate system with origin located at one of 
the points, for instance, the point Po. 

The equations of motion in this case are — by 


=1 
ats nis Sf . kt Ave. Seon misty 
j=) j=) ; 
Ix 
n—l H—] 
” Sk 16 my (% — %) a —%) me 
fo CS 1.11 
rr 
a) t a—l 
C == m 
ad = — mM, ~— + But Aaa 7D as mie 
Ist 


{=1, 2, ..., n~I 
where &,, 7,;, ¢, are the coordinates of P;. 


We shall also find it convenient to write these equations 
in vector form. Let i, j, k be unit vectors directed along 
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the &, , ¢ axes, respectively, and r;; the vector joining 
P;toP;. Then multiplying equations (1.11) respectively 
by i, j, and k and adding, we obtain 


a—) 
fe =m 84 Sim tS Sm, 2 
dt er) Fan _ ro 


fan, — os w=] 


We have just considered the motion of the system of 
particles relative to Py. However, we could consider the 
motion of all the particles relative to any particle P,. We 
may therefore write the following equations: 


a—1) ee | 


Pry i Eh s my tS Fyk 
act rot a : r »: me ys 
ke j= ‘J j=l J 
wih y (1.12) 
4=0, 1, eoop n—) 
k=0, 1, eaep ry | 
ighk 


4. Let us now consider the moment of inertia of 2 
gravitating particles: 


ae? a (1.13) 


(=1 


Differentiating this twice with respect to the time ¢ and 
taking into account that 


R= PRTH 


we obtain 
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ey" = 2 Sin, (x0! + yy! +22) +5; Pye +2 *). 


ar2 
éanl 


By making use of (1.2) and (1.4), we can write this equa- 
tion as follows: 


a7 =23 (a dey tM [ an t 9a; gz) +40 +2). 


Now, the potential function U is a homogeneous function 
of the coordinates of degree —1; therefore 


(a2 * ox, ae, + Y Wi 5 +4G)=— 


éun] 


and hence, 


a = 2(U + 2H). (1.14) 


This equation is called the Lagrange-Jacobi equation and 
plays an important role in celestial mechanics. It was 
obtained by Lagrange in the case of three bodies and was 
generalized to 7 bodies by Jacobi. 


5. We next derive an auxiliary equation which allows 
us to express the moment of inertia in terms of the dis- 
tances between the bodies. 

Consider the distances of the bodies Pi, Ps,..., P, 
from the center of mass of the system and the distances 
between the bodies. The square of the distance of P; from 
the center of mass is 
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ri =xit+yi tz, (1.15) 
and the square of the distance between P; and P; is 
735 =(%3—% 3)? +(91 95)? +(2;—2;)?. (1.16) 


We previously used }°;; to stand for the summation over 
all distinct pairs 27 (¢ + j and # not considered to be 
distinct from 2j); we now let 3-*;; denote the summation 
over all possible pairs 77 (any 2 is taken in combination 
with any 7). 

We then have the following obvious equations: 


S mm, (2, — 2,)* = > mM, (44— xP = 
ry) 6j 


oe R 
| 
—o_ 2 mm; (z,—2z,? = 
ém1 J=1 
i] ® 
4 >> 
— my > > mm jx + > m,m x4 — s Sires 
s=1 j=l =l j=l ‘=1 J=1 
M a 
=z Mt mjx's — » mH, S mM yX4, 
s=1 j=! s==1 I=! 
where . 
M=>m, 
=} 


However, since 


Zmat= 3 mz, 
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and 


by the linear momentum integrals, we find from the above 
equations that 


i] 
M > mz; = 2, mm, (%,— 2,)*. (1.17a) 
(=] 4J 
In an analogous way, we obtain the formulas 


" 
M> myi _ >, mm, (Ye yy)’ (1.17b) 
(=I $3 


M >, mzi = DY man, (2, — 2,) (1.17c) 
(=) ‘J 


Adding (1.17a), (1.17b) and (1.17c) and taking into ac- 
count (1.15) and (1.16), we find 


M Langa Zmmyiy (1.18) 
Let us introduce the notation 
Po pa m,m is. 
Then by (1.13) and (1.18), we have 


__ 1 7s 
L=ael. 


Thus, the Lagrange-Jacobi equation is expressible in the 
form 
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5 oF =2(U +24). (1.19) 


6. The energy integral and Lagrange-Jacobi equation 
permit us to establish certain properties of the distances 
r;; as time. becomes negatively or positively infinite, 
without the use of additional analytical tools, for the 
general case of a system of ” gravitating bodies. It is true 
that in this way we obtain only the simplest theorems on 
the final motion which already were given by Jacobi. 
However, the great importance of these theorems is that 
they reveal the proper way to formulate the problem of 
final motions, as well as indicate what general course a 
more detailed investigation should take. 

We first consider the properties of the distances 7;; 
in a system of ” gravitating bodies with a negative total 
energy. 

THEOREM 1.1. If the constant H in the law of conservation 
of energy is negative, and if as [> (t»—~) no collisions 
occur between the bodies, then for all t>0 (t<O) the function 


re(t)= min {r; 3} 
tj 
possesses the following properties: 
(1) For all t>0(¢<0) 
M3 
rea 


(2) For any T>0 no matter how large and any €>0 no 
matter how small, there exists a time t’ with |t/'|\2>T such 
that 
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M2 
relma: 
where : 
f= > Mm, 
(=! 


Let us prove the first statement. Since the kinetic 
energy of the system cannot be negative, it follows from 
the conservation law that 


U+H>0 
or 


Sy 2H 


rey 


We can strengthen this inequality and write 


Aw mm aH (1.20) 


However, it is not difficult to see that 


= M2 
> mmy Q-z 
tJ 


Thus without weakening inequality (1.20), we may write 


ra > 2 |AI, 


and hence 


MQS APE : (1.21) 
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We now prove the second statement. We show that for 
any positive 7, no matter how large, we can find a ’? >TZ 
for which 


U+2H>—e. 
Assume the contrary. Then for some value of T we have 


U+211< —e (1.22) 


for all values of t>T. 
Integrating the Lagrange-Jacobi equation (1.19) twice 
with respect to time from T up to ¢>T7, we obtain 


tt 
2(t)= 12 (T) + 22'(T)(t—T) + 2M [[(7 +24) dt dt. 
rT 
and hence by (1.22) 
Pt) <12(T) +I" (T)(t —T) — Me (t — TY? 


We see from this that [?(t)->—© ast»-+o. Therefore 
there must exist a time t>0 such that J?(7)=0, i.e., a 
collision between the bodies occurs. This contradicts the 
condition of the theorem. 

Thus, there exists a t/>TZ for which 


U(t')+2H>—e, 


U(t')>2|H|—e. 
Strengthening this inequality, we have 


M2 
2r,(t’) > | =e, 
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or 
M2 
/ ee erg 
ret) < Gyre ° 
and so for ¢>0 the theorem is proved. For ¢<0, the 
theorem is proved in a similar way. 

The fact that when JZJ<0, some of the distances can 
increase indefinitely with the time, does not contradict the 
theorem just proved. However, it follows from the theorem 
that this cannot happen to all the distances. Thus when- 
ever H<0, it is impossible for a system of gravitating 
bodies to be completely scattered either as t>o or as 
[-->— co. 

We now consider a system of gravitating bodies for 
which the total energy is positive, i.e., H>0. 

THEOREM 1.2. (Jacobi’s Theorem). If the constant H in 
the energy tntegral is positive, and if as t->~(t>+—@) no 
collisions occur between the bodies, then at least one of the 
distances 1; ; increases indefinitely. 

Proof: We start with the Lagrange-Jacobi equation 
in the form (1.19) and integrate it twice with respect to 
time from ¢=0 to +>0. We thus obtain 


tt 
P() = B+ t+ aM [[(U + 2H) deat, 
nn 


where J% and J’ are the respective values of J? and Jz 
at t=0. Now, since U>0, it follows that 


te 
P(t) > P+ +4M | | Hae dt, 


or 


P(t)> [24 [t+ 2MHe. 
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When H> 0, we find from this inequality that 
PQ) > asto>+o (t+—~), 


This is only possible if at least one of the distances s, , 
increases indefinitely as i>-+ © (t>—«). The theorem 
is proved. 

This theorem is not incompatible with the fact that, 
when H> 0, instances are possible in which all of the dis- 
tances between the bodies increase indefinitely as to, 
or as +>—o, or else in both cases. However, the 
theorem is also not incompatible with the fact that, when 
H>O, there may be subsystems in which the distances 
between the bodies remain bounded at >, orasi>— , 
or as approaches both +. Finally, the fact that there 
may be complete scattering of the system as the time 
increases in one direction coupled with incomplete scatter- 
ing as the time increases in the opposite direction is not 
incompatible with the theorem. 

In conclusion, by comparing Theorems 1.1 and 1.2, 
we can state that, when H<0, the energy in the system is 
insufficient for it to be completely scattered, and when 
H>0, the energy is so high, that at least a partial scatter- 
ing of the system is inevitable. 


CHAPTER 2 
Method of Dimensional Analysis 


1. The method of dimensional analysis, which has 
lately been applied rather widely in diversified problems 
of mechanics and physics, can also be used in the qualitative 
analysis of the motion of systems of gravitating bodies. 
However, it must be borne in mind that dimensional 
analysis is incapable of yielding complete enough results 
for the -body problem, though in certain instances it 
enables important conclusions to be obtained. 

Our purpose is to consider means of applying the meth- 
ods of dimensional analysis to find criteria for the occur- 
rence of various types of motions in a system of 7 grav- 
itating bodies. In those instances where dimensional 
analysis enables us to obtain such a criterion, we find that 
it is expressible exactly to within an unknown dimension- 
less factor. Thus, the methods of dimensional analysis 
do not lead to a complete criterion, but only give an 
indication of its analytical form. Nevertheless, such re- 
stricted results are of importance since they yield pre- 
liminary knowledge and give us a correct orientation at 
the initial stages of our investigation. 


2. That dimensional analysis is capable of indicating 
the analytical form of the dependence of two dimensional 
quantities springs from the fact that objective physical 
relationships cannot depend on the choice of the system 
of units. However, such relationships between quantities 
are not arbitrary and are determined by the dimensions 
of these quantities. 

Dimensional analysis is usually applied to some rela- 
tion expressible in the form of an equation or identity. 

15 
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What is somewhat new in the application of the method 
of dimensional analysis to the qualitative study of the 
n-body problem is the necessity for considering inequalities. 
For our further purposes, the following reasoning will 
suffice. 


(1) Let SW be a set of systems.of basic independent units 
Gi, Ja) - + +> Qn 
such that a transformation to any other system of units 
911 G21 +++ Un» 


also belonging to 9, is effected by means of the formulas 
= Odi (2.1) 


where a; are arbitrary positive numbers. 
(2) There exist arbitrary positive quantities 


Qo, 1, Qz, soy Ons 


which are homogeneous in the sense that the values of these 
quantities: 


Qo, O'1, Q'2,--- 5 Os 


after transforming the basic units according to (2.1), are 
related to their previous values by the equations 


O';= afsta$i2... atinQ;. (2.2) 


The exponents a;, are the powers of the dimensions; if 
they are all zero, then the Q; are dimensionless. It is clear 
that the product of homogeneous quantities is again 
homogeneous. 
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(3) The quantity Qo is related to Qi, Qo, ...,Q, by 
the inequality 


Qs fQu Qe, oes Qn). (2.3) 


(4) The relation Qo<f(Q1,Qo, ..., Qy) does not depend 
on the units of measurement in the sense it is satisfied after 
every transformation of the basic units by means of (2.1). 

(5) All the basic units enter into the expressions for 
the quantities Qo, Qi, ..-, Q,. 

THEOREM 2.1. Jf sanasions (1) through (5) are satisfied, 
then we have 


Qu $ CQOTOF . .. Qin, (2.4) 


where the dimensionless constant C=f(1,1,...,1), and the 
exponents x; are uniquely determined by condition (4). 
Proof: We take the quantities Q;, Q2,..., Q, and 
transform to new units; we choose the scales a1,a@o, .- -) Qs 
so that Q}, 02, ..., Qf each has the value 1 in the new 
system of units. Imposing this condition on (2.2), we 
obtain 
afilags2... Gene ; (2.5) 
j 
Solving these equations, we find that the a, are given by 


1 

~ QiQhie...QBin’ (2.6) 
where b;, are constants whose values we shall not deter- 
mine. (We might suggest that the reader solve equations 
(2.5) for the simplest case of n=2 if he is interested in the 
details of the derivation of (2.6)). 

Let us now write (2.2) for the particular case where 
j=0: 
Qo = ara .. . amo: 
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Substituting in this the value of a, from (2.6), we obtain 


»_——__Qo 
CO Of Oia. Oe’ me 
where for the moment we do not determine the exponents 
Nyy Lay - oo y Ly. 

Inequality (2.3) does not depend on the units of 
measurement, and it is therefore also satisfied after going 
over to the units gi, q2,.-., 9%; thus 


QoSf(1,1,..., 1). 
Setting 


C=f(1,1,..., 1) 


and substituting in the last inequality the expression for 
Qo given by (2.7), we obtain 


___Qo _¢¢ 
OF10%2...Q%n ~~? 
and hence 


Oo < CQ%10%2 oe - OFn : 


Inequality (2.4) has thus been obtained. Now, according 
to condition (4), the x; are such that this inequality is 
dimensionally homogeneous, i.e., the powers of the dimen- 
sions on the right and left-hand sides must be the sarre for 
each basic unit. Thus equating them, we obtain equations 
from which the x; may be determined. 


The theorem is proved. 

3. Let us consider a system of m gravitating bodies 
Po, P1,..., Px—1 with masses mo, 11, ..., M1. Let 1; ;(t) 
denote the distance from P; to P; and 1; ;(t) the derivative 
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of this quantity with respect to time. We shall call the 
latter the radial velocity. 
Let us introduce the notation 


p(t)= min {7r;;(t)}, 
o(t)=min {rj ;(t)}. 


We shall assume that at the initial time ¢=0 the condition 


o(0)>0. (2.8) 
is satisfied. 

Under this assumption, there are two possibilities: 

(a) it is possible to give a positive number 7 such that, 
for all ¢ in the half-open interval [0,7), the distances 
between the bodies increase and from the moment of time 
T at least one of the distances ceases to increase; 

(b) or if o (0) is so large, and therefore also the initial 
values 7;;(0) of all the radial velocities, that the mutual 
attraction of the bodies of the system is not great enough 
to contain the system, then all the 7;,;(¢) will approach 
infinity as t->o. 

However, for exactly what value of o (0) will this 
phenomenon occur? This, of course, will depend on the 
initial state of the system, on the parameters of the system 
and on the parameters which characterize the phenomenon 
of gravitation in general. 

Let us try to clarify exactly which parameters are of 
great importance in the problem under consideration. 

From general physical notions about the force of 
gravity, it follows that the smaller the value of p(0), i.e., 
the closer the bodies are at the initial time, the greater 
must be the value of a (0); therefore, there is no doubt 
that one of the parameters is p(0). 


20 QUALITATIVE METHODS IN THE MANY BODY PROBLEM 


The masses 219,71, ... , M,—; are also of considerable 
importance. For, the larger the masses of Po, Pi,...,P,—1, 
the stronger will be their mutual attraction and, with all 
other conditions being equal, the greater must be the value 
of 7 (0). Thus, taking into account the masses will yield 
more parameters. However, we can avoid the large number 
of parameters entailed by the masses if we remember 
that we are going to obtain sufficient conditions in the 
form of certain inequalities. Therefore for the purposes 
of simplifying the problem, we can use an overestimate of 
the role of the masses. As a parameter characterizing 
mass effects, we take 


b= (mo mM, ..., My—-1), 


which is some function of the masses having the dimension 
of mass. 


Finally, gravity is of essential importance, and there- 
fore the gravitational constant yy must also be one of the 
parameters of our problem. 


The parameters y, », p(0) and oa (0) constitute a 
complete set of parameters for the phenomenon under con- 
sideration. In fact, for given values of p(0) and yz, if 
o (0) is sufficiently great, then the bodies of the system 
will all be scattered, regardless of what values the other 
unaccounted for parameters may have. “Thus we con- 
clude that the required criterion must have the following 
general form: 


a(0)>f(y, wu, p(0)), (2.9) 


where the function f is momentarily unknown. We must 

now determine the functional form of f(7,“,p(0)). 
Consider arbitrary but fixed values of the parameters 

7, wu, and p(0); leta.(0) denote the greatest lower bound 
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of the set of values of ¢ (0) for which the system of bodies 
is completely scattered. This lower bound obviously exists 
because of inequality (2.8). It is further obvious that 
the inequality 


o%(0)< f(y, w, p (0)) (2.10) 


is satisfied for all values of y, », and p(O). 

Let us take the system of units ordinarily used in 
problems in mechanics in which the basic units are time 7, 
length Z, and mass M. The parameters in which we are 
interested have the following dimensions: 


[p(O)J =L, 
_— = (2.11) 
[yJ=LT2?M—. 


We thus have four parameters, three basic units, and 
the inequality (2.10), i.e., our conditions correspond to 
those of Theorem 2.1, and to apply it, we must set n= 3. 
By this theorem, we then have 


F4(0) S Cy*n* p(0)*, (2.12) 


where the exponents x, y, and gz are for the moment un- 
known. To determine them, we equate the dimensions 
of the left-hand and right-hand sides of this inequality. 
This yields 
LT = T—**L3* M-* M°L', 
or 
LTA = L3*+27—-2* Yyr-=., 
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Equating exponents of like dimensions, we obtain the 
equations 


2x=1, 
y—x=0, 
3x+2=1, 
whose solution is 
| i \ 
i=7) y=’ Jee Ti 


Substituting these values of x, y, and z in inequality 
(2.12), we obtain 


%(O<C vine (2.13) 


However, if our system of units is the one usually used in 
theoretical investigations and in which y=1, then the 
criterion becomes 


1 (0)<C VV OR (2.13a) 


Our above result implies that 
ee Biba 
fy, a p(0))=C > (0) ” 


and thus (2.9) becomes 


)>CV 45 (=): (2.14) 
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The value of C and the form of the function expressing 
M in terms of m,7m, ..., and m,_, cannot be determined 
by means of dimensional analysis. Thus dimensional an- 
alysis does not yield definitive sufficient conditons for the 
occurrence of the type of motions of n gravitating bodies 
in which we are interested; it merely determines the 
analytical form of the criterion. Later on, we shall see 
that the results obtained are confirmed by the direct 
analysis of the equations of motion. 

In conclusion, let us examine the methodological as- 
pects of our above investigations. 

We obtained inequality (2.9) on the basis of general 
physical considerations. This part of the investigation 
does not concern dimensional analysis. Figuratively 
speaking, inequality (2.9) is not within the ‘‘competence”’ 
of dimensional analysis, and it cannot be “‘responsible’”’ 
for this inequality. What merely follows from dimen- 
sional analysis is that if (2.9) is in principle correct, then 
it has the form (2.14). The question as to whether or 
not inequality (2.9) is the correct form of a sufficient 
condition for the occurrence of different types of motion 
cannot in principle be answered by dimensional analysis. 

In this connection, it is interesting to note the follow- 
ing. Our reasoning cannot be applied to the case where 
7 (0)<0. In fact, if we wrote the inequality (2.9) in 
the form 


lo (O)|>f(Y, u, p(0)), (2.15) 


we could, at first glance, repeat our argument and obtain 
a similar but no more general result. However, in reality, 
this is not so. For, in assuming that o (0)<0, we would 
not have the right to suppose that the problem is solved 
by inequalities (2.9) and (2.14). When o(0)<0, it is 
possible for the gravitating bodies to approach one another 
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very closely, and the resulting gravitational interactions 
could then lead to such an exchange of energy between the 
bodies that a portion of them could go off to infinity with 
increased velocity, whereas the energy of the remaining 
bodies would be insufficient for them to be scattered. In 
this case, it is not sufficient to consider merely the 
parameters ¢ (0), e(0), u and y in formulating conditions 
for the scattering of gravitating bodies, and we require 
a more detailed set of parameters that characterizes more 
fully the initial state of the system of gravitating bodies. 
Therefore, when ¢ (0)<0, the relation (2.14) must either 
be replaced by one containing a fuller set of parameters, 
or it must be supplemented by relations containing other 
necessary parameters, or else our considerations must be 
limited to such a set of initial states for which the un- 
desirable close approaches cannot occur. 


Let v;; be the magnitude of the velocity in the relative 
motion of P,; and P;, and let 


v(t) = min{y, ;(£)}. 


It is clear that if motions of gravitating bodies are to 
occur for which p(t)—>© as >, it is not at all necessary 
for the radial velocities to be sufficiently great. Such 
motions can occur if the absolute velocities in the relative 
motion of the bodies are sufficiently great. Nevertheless, 
we cannot substitute v(0) for o (0) in inequality (2.9), 
even though the dimensions of these quantities are identi- 
cal. This is due to the fact v(0) is always positive, whcreas 
it is possible for v(0)=|oa(0)| when ¢ 0)<0. But we saw 
above that a criterion of the form (2.9) does not apply to 
this particular case. Therefore, any criterion for the 
occurrence of motions in which 7; ;(t)—>© asi» cannot 
contain merely (0), (0), uw, and +, but must contain 
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still other parameters that, to one extent or another, 
explicitly or implicitly account for the direction of the 
velocity. As we shall see later on, this may be achieved 
by using v(0) and o (0) together since their difference, to 
a certain extent, characterizes the direction of the velocity. 
In this way, one may obtain a criterion which is applicable 
for a sufficiently wide range of values of ¢ (0)<0. 

We have limited ourselves to a comparatively simple 
application of the method of dimensional analysis. How- 
ever, this application makes rather well-apparent its 
characteristic features and the possibilities of its use. 

In summing up, we must admit that the method of 
dimensional analysis requires great care and thought, as 
well as a deep and penetrating intuition into the essence 
of a problem. Nevertheless, dimensional analysis can 
be earnestly recommended if it is simply regarded as an 
exploratory or research means in theoretical considerations. 


CHAPTER 3 
Method of Continuous Induction 


1. We shall consider aset SN whose elements are time- 
dependent processes or mechanical or physical systems S 
whose states vary with the time ¢. 

Let I’ be some property of S. In what follows, we 
shall make use of the following definitions. 

(1) The property IT issaid to be tnductive on the set NM 
if for every SCO we have the following : from the assump- 
tion that the property I’ is satisfied for all values of ¢ 


such that 
0<i<r, 


it follows that it issatisfed for =r no matter what 7 is. 

(2) The property I issaid to be continuous if from the 
condition that it is satisfied ata certain moment of time 
ty, it follows that it is satisfied in a certain neighborhood 
of 4, i.e., itis possible to imda number 5>0 such that [ 
is satisfied for all such that 


[i—tl<d. 


THEOREM 3.1. (On Continuous Induction). Let the 
property T be continuous and inductive onthesetm. If T 
ts satisfied at the initial teme t=0, then tt 1s satisfied for 
all values of t>0. 

Proof: Wecall the time ¢ regular if for all ¢ such that 


0 <0’ <8, 


the system S possesses the property I. Let Q denote 
the set of all regular times. From the continuity of [I 
and the assumption that I is satisfied at ‘=0, it follows 
that the set Q is not ernpty. 

27 
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We prove the theorem by contradiction, and we assume 
that I is not satisfied for all values of t20. Therefore, 
Q is bounded from above; let r=sup Q. Since 7 is the 
least upper bound, it follows that 7 is satisfied for all 
values of ¢ in the interval 0<i<7. But then since I 
is inductive, it is satished at time 7. Furthermore, the 
property I’ is continuous. Therefore, we can find a 
5 >0 such that I is satisfied for all ¢ for which 


|i—T |< 6, 
But then I. must be satisfied for values of ¢ for which 


O<tgr+ 2 
2 
Hence, it follows that the time 
i) 
T1=T+ 5 >T 


is regular. However, this contradicts the fact that 7 is 
the least upper bound of the set Q. The theorem is proved. 
2. To apply the Theorem on Continuous Induction, 
we must have criteria for when some particular property 
is actually inductive. We shall now show that such 
criteria can be given in certain sufficiently general cases. 
Let 
W1, WO, 20, Wy 


be some set of parameters of the system S depending on 
the time which, while not completely determining the 
state of the system, characterizes it in certain respects. 
We denote the set of functions 


G1, Wa, .. «, W, 


by {w}. Moreover,let the system depend on » parameters 
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Gy, Ag, ..-, Quy 


which are independent of the time, and which we denote 
collectively by {a}. 

Let A({a},{w}) denote an operator which, according 
to certain rules, sets in correspondence to each set of 
values of the parameters dy, d2,..., @, and to each set of 
y functions w;, wz, .., w, a continuous function W(t) 
defined for all #20. 

We first consider the case where 1, we, ...,w, are 
unknown functions, and our knowledge. of the system S 
is restricted to the information that these functions and 
the parameters aj, dz, . . ., a, satisfy relations of the form 


w;(t)2A;({a}, {w}). (3.1) 


$m1,2,..;, y 


Here, A; are operators which are defined by means of 
their analytical expressions. 
Let 


Filt), felt), see St) 


be suitably-chosen continuous functions of t defined for 
all 20 which we shall call comparison functions. The 
property I" will be that the following inequalities are 
satisfied : 

wi) >Si(t). 

om. 2. wey ® 

THEOREM 3.2 (First Theorem on the Inductiveness of 

Inequalities). If a system S sattsftes the operator inequalities 
(3.1), and af for any t20, from the assumption that the 
inequalities 

w(t) >f;(t) (3.2) 

é—1, 2, y 
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We prove the theorem by contradiction, and we assume 
that T' is not satisfied for all values of 20. Therefore, 
Q is bounded from above; let r=sup Q. Since 7 is the 
least upper bound, it follows that 7 is satisfied for all 
values of ¢ in the interval O<i<r. But then since I° 
is inductive, it is satisfied at time +r. Furthermore, the 
property I is continuous. Therefore, we can find a 
6 >0 such that I is satisfied for all ¢ for which 


jti—Tr |< 6, 
But then I’ must be satisfied for values of t for which 


O<igrt ud 
2 
Hence, it follows that the time 
5 
T1=T+ Zi >T 


is regular. However, this contradicts the fact that 7 Is 
the least upper bound of the set Q. The theorem is proved. 
2. To apply the Theorem on Continuous Induction, 
we must have criteria for when some particular property 
is actually inductive. We shall now show that such 
criteria can be given in certain sufficiently general cases. 
Let 
W1, Wa, 220, W, 


be some set of parameters of the system S depending on 
the time which, while not completely determining the 
state of the system, characterizes it in certain respects. 
We denote the set of functions 


W1, Wa, ..-, W, 


by {w}. Moreover,let the system depend on » parameters 
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G1, Go, .--, Us 


which are independent of the time, and which we denote 
collectively by {a}. 

Let A({a},{w}) denote an operator which, according 
to certain rules, sets in correspondence to each set of 
values of the parameters aj, @,. . ., @, and to each set of 
y functions w, we, .., w, a continuous function Y(t) 
defined for all #20. 

We first consider the case where w:1, we, ...,w, are 
unknown functions, and our knowledge. of the system S 
is restricted to the information that these functions and 
the parameters aj, a2, . . ., a, satisfy relations of the form 


wi O24; i Ug}, {o}). (3.1) 


Here, A; are operators which are defined by means of 
their analytical expressions. 
Let 


fit), Salt), ss ft) 


be suitably-chosen continuous functions of t defined for 
all £20 which we shall call comparison functions. The 
property I" will be that the following inequalities are 
satisfied : 


THEOREM 3.2 (First Theorem on the Inductiveness of 
Inequalities). If a system S satisfies the operator inequalities 
(3.1), and af for any 720, from the assumption that the 
inequalities 

Oj (2) > Fit) (3.2) 
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are saltsfied in the half-closed interval O<t<7, tt follows 
that the inequalities 


A ;({a}), {o})2A:Ca}, {f}), (3.3) 


are satisfied, then (3.2) 1s inductive on the set St of systemts 
S for which the inequalities 


i({a}, Fil) > (3.4) 


a=1, 


are satisfied for all t20. 
Proof: Let us write inequality (3.1) in the form 


w()—f(t)24; ({a}, i}) —f,(t). (3.5) 


$™1, 2, 206, 
Choosing an arbitrary T >0, we assume that 


w (t)>f;(2) for ak (3.6) 


swt, 2, 20, 


Then considering (3.5) for values of ¢ in the interval 
0<i<7 and using (3.3), which holds because of (3.6), we 
can write 


0()— fil) 24 s(ia}, {f)—fi(d. 


We now let ¢—> 7 taking into account (3.4) and the fact 
that w;(¢), f;(¢), and A;({e}, {w}) are continuous func- 
tions of ¢, and we obtain in the limit 


Wy (7) Flr) 2-4 Ua}, )—-fi(7)>0, 


or 
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w3(7T)>f;(r). (3.7) 


> ay eee 


We have thus shown that the inequality (3.7) follows 
from (3.6), and this proves the theorem. 


3. We now consider a more complicated case of the 
Theorem on the Inductiveness of Inequalities which we 
shall prove in two parallel versions. 

Suppose we have two sets of parameters of the system 
S which depend on time: 


Wi, Wa, -.-, Wy, 


61, 50, ..., 5, 


which, though not defining the state of the system com- 
pletely, characterize it in certain respects. As before, 
suppose the state of the system also depends on » para- 
meters a1, G2, . . ., @, which are independent of time. We 
consider the case in which w, wo, ...,w,and 61, de,..., 5, 
are unknown functions, and our knowledge of the system 
S is restricted to the fact that these functions and the 
parameters a1, de, . . ., a, satisfy relations of the form 


5 (t)<A,({a}, {o}) 
(8;(2)2A,({a}, {w})). 


tel. 2). occa: ¥ 


(3.8) 


Here, A ; are operators defined by their analytic expressions. 


Let 


v(t), 2 (2), veey y(t), 
fi(2), fr, a f(t) 


be suitably chosen continuous functions of the time # 
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which are defined on the range ¢>0, and which we call 
comparison functions. We shall assume that they pos- 
sess the following property: From the assumption that 
the inequalities 


G, (t) < so (t) 
(3,() >e(t), (3.9) 


4=1, 3, eoce ¥ 


are satisfied in the interval 0<i<7, it follows that the 
inequalities 


u(t) Df, (2). (3.10) 


é=1, 2, eee ¥ 


are satisfied for all values of t. 

Although we do not know the functions w;(¢) and 
5;(¢), we can, however, often establish such relations 
from the physical or mechanical meaning of these func- 
tions or from other considerations not requiring the know- 
ledge of the analytical form of these functions. 

We shall take I’ to be the property that the inequali- 
ties (3.9) are satisfied. | 

THEOREM 3.3 (Second Theorem on the Inductiveness 
of Inequalities). Jf the system S satisfies the operator 
inequalities (3.8), and tf for any Tt >0, from the assumption 
that the inequaltties (3.10) are satisfied in the interval 
O<t<r tt follows that the inequalities 


A,({a}, {o))< Ay ({a}, {(f}) 
(A,({a}, (o}) > 4,({a}, (1)))> 


‘=1, 3,. ood ¥ 


(3.11) 
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are satisfied, then property (3.9) (and (3.10) which follows 
from it) ts inductive on the set M of systems S for which the 
tnequalities 


A, ({a}, (/}) — 4 (¢) <0 
(A, ({a}, (}) — 9, (t) > 0). 72) 


‘=l, 2, eee 


are satisfied for all t20. 
Proof: Let us write inequality (3.8) in the form 
8, (2) — 9 (4) < A, ({a}, {w}) — % (¢) 
(3,(t) — 9 () > A,({a}, {o})—9eft})  -13) 
Cee Be er 


Choosing an arbitrary tT >0, we assume that 


95 (4) << 9 (2), 
O0<t<c: 


(3,(t) > 9 (2)). 


é=21,3, ..03 


(3.14) 


Then for these same values of ¢, the inequalities 


wo, (f) > Ie (?). 


(= 1,2, Cr 


(3.15) 


will be satisfied. Considering (3.13) for values of ¢ in the 
interval O<i<7r and making use of (3.11), which hold 
because of (3.15), we can write 
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We now let t->7 taking into account (3.12) and the fact 
that w ;(t), 5 ; (2), f:(), y;(t), and A ;({a},{w}) are con- 
tinuous functions of /, and we obtain in the limit 


(— OSA ({4}, (1}) — (7) CO 
(5 ()— el) > A ({a}, (1) — 94 (2) > 0). 
6mm], 2, os os 


or 


8, (t) << 9, (t) 
(3, (<)> 9, ()), (3.16) 


‘=1,2, ....% 
From this, in turn, it follows that 


w, (%) > fr (2). 


é=}, 2, eee *¥ 


We have thus shown that inequality (3.16) follows from 
(3.14) with which the theorem is proved. 

4. We now consider some applications of the Theorems 
on Continuous Induction and the Inductiveness of In- 
equalities to the qualitative analysis of gravitating bodies. 

We begin with the first Theorem on the Inductiveness 
of Inequalities. 

Let us take a system of n gravitating bodies 
Po, Pi, .. ., Py. with masses mo, m, . . ., M,—1, and let 


OI gp my 


We shall make use of Jacobi coordinates. To begin with, 
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we consider the motion of P; relative to Py. Differen- 
tiating the identity 


ae =—3 1 a ue oad Ss 
twice with respect to ¢, we obtain 


; if 2 12 
rola to = bey amber yc’, 
or 
" nn 1" Se ” 
"oo = ae a ro hy oi + 
are 2 2 "9 
Goth th ly 
_ 


ro. 


Taking into account that 


12 2 2 7) 
Bt +o) — 254 >0 
Fo. ot 


we find that 


r ge 


Patt Poi 


n+ LG. 


Fo. 


Using the first of equations (1.7), i.e., for 7=1, we obtain 


” 1 aU 
ra? etc ae et Ti: %) 


Fou 


36 QUALITATIVE METHODS IN THE MANY BODY PRORLEM 


Now, let / denote the direction of the line joining Po to 
P,. Then the last inequality may be written in the 
following form 


” { aU 
ro (3.17) 


Differentiating the potential function U in the direction 
l using (1.1), we have 


au S mymy Orsy 
on 
ed 


From simple geometrical considerations, it follows that 


Ores 
els 
and, therefore, 
au miyIny 
~2-y——- (3.18) 


Inequalities (3.17) and (3.18) imply that 


. | 1 myn; 
Ho] om P; R 
J J 


We have obtained this inequality by considering how 
P; moves relative to Py. However, we could have con- 
sidered the motion of any point P, relative to any other 
point P,. Thus, we have 


1 WV mimy 


“ee aaenal ry 
‘> 


ww” 
a 


Finally, integrating this inequality twice with respect 
to time from ¢=0 to ¢>0, we obtain the following system 
of integral inequalities: 
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st 
- mem 
rey (t) > rr, (0) +7}, (0) ¢ — nog j | = dt dt. (3.19) 
00 
We now use these inequalitics to find sufficient con- 
ditions which ensure that the distance of each body from 
the remaining ones will become infinite as >, in other 
words, that the system is completely scattered. Such 
will be the case, at any rate, if the inequalities 


Pye (0) rn) 
I(t) > FE : 


are satisfied for all #20. 


To find the conditions, it is most convenient to use 
the first Theorem on the Inductiveness of Inequalities 
(Theorem 3.2). 

For the functions w;(t), we take the distances r,, (t) 
and for the parameters a,, the initial values of these dis- 
tances and their derivatives, 7,,(0) and r’,,(0), and the 
Masses 770, 721, . - +» Myz—1. 

In accordance with inequality (*), it is natural to 
choose the comparison functions to be 


a (U) ry (Y) 
sae 


We define the operators A,, using (3.19) as follows: 


hea (j= — 


=p, (0)-++ 7, (0) ¢— isan) 

i) 
—_— f'dt dt. 

dais ‘J ee ry 
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Finally we take I’ to be the property that the following 
inequalities are satisfied: 

tas(t)>fi s(t). (3.21 ) 


The conditions of Theorem 3.2 are satisfied since, from the 
assumption that (3.21) are satisfied in the interval0<i<T, 
it follows that 


Therefore, by Theorem 3.2, T is an inductive property 
for those sytems of gravitating bodies for which the 
inequalities 


1 (0) uP 


+ —— 
$ ft 


1 mymy 
= GE ee eee ee gee ee dt dt 0. 3.22 
Hee > | (4 4 rg) ‘) pone 
0 «0 2 2 


are satisfied for all ¢20. These inequalities should be 
considered as defining the appropriate values of 7; ;(0) 
and r; j (0). 

Evaluating the integral on the left-hand side of (3.22) 
we find 
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tet 

\\ ay ta tt 
ri, ( Pej 

Jo) ( a 


4mymy ; 4mm, In (1+ r;; (0) | 


rey (0) yy (0) ry (0) ry (0) 


—e 


Substituting the value of the integral in (3.22), we obtain 
after simplification 


| ry, (0) 1 4mm, 
| 2 Mee > ry (0) 7, 5 5 


é 
(3.23) 


ee Py, (0) ae a — >> isi Ee in( 1+ ray (0) (0) eet) > 0. 


Mes “ss ry | (0) Peg (0) 


The second and third terms on the left-hand side are 
positive for every value of t>0, provided that 7’; ;(0)>0. 
Thus, (3.23) will be satisfied for all positive ¢ if the co- 
efficient of ¢ in the first term on the left-hand side is 
positive, i.e., if 


Pe, (0) _ ft >) 4mm, 
2 Pee “G 7449) ray (0) 


When these conditions are satisfied, inequality (3.21) 
will hold for all positive t, provided (3.21) holds for :=0. 
But then (3.21) is satisfied at t=0 since 
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fry (0) = THO) 


We can thus formulate the following theorem. 
THEOREM 3.4. If the conditions 


r (0) >, 


rg) (0) ~—>, 4m, sm y >0 (3.24) 
Set r, 5 (0) 70 


are satisfied at the inittal time t=0 for a system of gravitating 
bodies, then the distances r;; between the bodwes tend to 
enfinity as [>o. 

5. The criterion (3.24) just obtained is rather simple, 
but, nevertheless, slightly cumbersome. However, by 
weakening it, we can reduce the criterion to a simpler 
and more transparent form. 

Let us introduce the notation 


p(= mn (743 (t)}s 


s drjy (t) 
off = min} 
( ) Y, dt 9 
ngIn 7 
M*—-2__ 
min wy © 


It is not difficult to see that conditions (3.24) hold if 


3(0)>0, 


at (0) > _ (3.25) 


Thus, if conditions (3.25) are satisfied, the distances 7; ; 
between the bodies tend to infinity as t>-~. 
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We might note that conditions (3.25) have exactly the 
same analytic form as those which were obtained in 
Chapter 2 using dimensional analysis. 

6. Let us now consider a system of gravitating bodies 
Po, Fi, ..., P,—1, the motion of which we shall describe 
using Jacobi coordinates. Set p=p,—1, where p,—.: de- 
notes the distance of P,,: from the center of mass of the 
remaining bodies. Choose a positive number a <1, and 
consider the sphere S(a) with center at the point P,—; 
and with radius R=ap. Thus, the center of the sphere 
moves together with the point P,_1, and its radius varies 
with time. We shall say that an a-approach occurs be- 
tween the body P,—: and some other body P; at time # 
if at this time P; lies inside or on the boundary of the 
sphere S(a). 

Differentiating the identity 


pr= FF tet o7 


twice with respect to the time, and carrying out the same 
kind of manipulations as was done for 7 in Section 4, 
we obtain 


ee gi (3.26) 


where / is the direction of the line joining the center of 
mass of the bodies Po, Pi, . . ., P,—2 to the point P,_4. 
Since in differentiation with respect to the direction / 
only the distances between P,_, and the remaining 
bodies vary, we have that 


n—3 
ba désen Or wd 


OL 


é=0 ri w= 


42 QUALITATIVE METHODS IN THE MANY BODY PROBLEM 


From simple geometrical considerations, it follows that 


or é,u—l 
| st | 1 
Therefore 
a= 
dU mgmy—} 
aan > x (3.27) 
{=O é, n—i 


If it is assumed that no a-approach between P,_, and 
the remaining bodies of the system occurs when /2 0, then 


lin > ap. 


Using this to strengthen the previous inequality, we obtain 
M 
pF, ’ (3.28 ) 
where 
a=) 
Me 2," 
Finally, if we integrate (3.28) twice with respect to time 


from t=0 to ¢>0, we obtain the integral inequality 


tt 
0) >eO)+e(O—Sl[F. 6.29) 
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THEOREM 3.5. If the tntital conditions 
p’ (0) >0, 


2 8M (3.30) 


hold at t=0 for a system of gravitating bodies Po, Pi, ..., Py—1 
and tf no a-approach of P,,—1 to the other bodies occurs for 
any value of t>0, then the distance of P,,-, from each of the 
other bodies in the system becomes infinite as ->o. 


Proof: The integral inequality (3.29) is of fhe same 
type as those in the system (3.19). Applying Theorem 
3.2 to it and repeating the argument that was used in the 
proof of Theorem 3.4, we obtain a proof of the stated 
theorem. 


7. We now consider some applications of the second 
Theorem on the Inductiveness of Inequalities. As a first 
example, it will be constructive to consider one of the 
cases already treated with the help of the first Theorem 
on the Inductiveness of Inequalities. Let us turn again 
to Theorem 3.5, and let us prove it in a slightly stronger 
form. In so doing, we obtain some results which we shall 
find useful below. 


THEOREM 3.5a. If the initial conditions 


p’(0)>0, 


‘ 4\f 
‘2 sates 


hold at t=0 for a system of gravitating bodies Po, Pi, ..., Px—1; 
and tf no a-approach of P,,—1 to the other bodtes occurs for 
any value of t20, then the inequalities 
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p> EO) | (3.31) 
p(t) >p(0) + Ors, (3.32) 


hold for all positive t, and, therefore, p(t)—>o as tro. 
This implies that the distance of P,-1 from each of the other 
bodzes in the system becomes infintie as t-> ©. 

To prove the theorem, we start with inequality (3.28). 
Integrating it with respect to time from #=0 to ¢>0, we 
obtain 


(2) >P()—S f 


Applying the parallel version of Theorem 3.3, we set 
8 {t) =p (?), 
w {t)=p (8), 
and we define the operator A by means of the equation 


t 
dt 

= MIF, 

For the comparison functions, we take 


_ ¢ (9) 
g(t)=>-, 


(0) 
f(t) =p (0) +=". 


It is first necessary to show that the conditions of 
Theorem 3.3 are satisfied for the case in question. With 


METHOD OF CONTINUOUS INDUCTION 45 


our choice of 5(¢), w(t), p(t), and f(t), it is obvious that 
(3.9) holds. The operator A also satisfies the required 
conditions since from the asumption that p(t)>/(¢) in 
the interval 0<t<7, it follows that the inequality 


d 


t t 
O- Blg>eO— ala 


is satisfied for these same values of ¢. Therefore, by 
Theorem 3.3, inequality (3.31), as well as inequality (3.32) 
which follows from it, are inductive properties for systems 
of gravitating bodies for which the inequality 


é 
ae (Q) Mee 
a) (+ (0) +e)’ = 


holds for all values of t>0. After integration, this last 
inequality becomes 


2 
“a Rey mw) 


at (» (0) + ~Z ) 
It is easy to see that when the conditions of the theorem 
are satisfied, this inequality holds for all values of #20. 
Therefore, under the conditions of the theorem, (3.31) 
and (3.32) are inductive properties, and since they hold 
at the initial time ‘=0, they are also valid for all #20. 
The theorem is proved. 
8. Let us now consider some further applications of 
the second Theorem on the Inductiveness of Inequalities. 
Suppose we have a system of gravitating bodies 
Po, Pi, ..., P,—1 with masses m, m1, ..., m,—1. To de- 
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scribe their motion, we make use of equations (1.12). 


Integrating these equations with respect to time from ¢=0 
to !>0, we obtain 


f 
ae | 
' r 
ry oo Oy (O) = -—m, 7 =e dt + a” fy dt -- 
Tok 
eur 
wd f 
| Fy9k 
= x | 73 dt. 
= jk 
j=l 
Ik a 


From this, it follows that 


t 
n=! 
; dt Sn dt 
nee |< J At j=90 Nt = 


jth 


ry (3.33) 
With a view towards applying Theorem 3.3, we set 
we (¢) = 7%, (t)s 
8 (= [re — FeO) 


and we define the operators A({a;,}),{w}) by means of 
the equations 


i n—1 ; n—l1 ; 
d dt de 
damm | = + Sm, | 4 >" age 88) 
" 


ss | an 
jp 
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We can thus consider (3.33) to be operator inequalities 
satisfied by the system of gravitating bodies. 
We now assume the initial conditions are such that 


Y (0) — 2ri, (0) > 0, 
(3.35) 


Uy (0) Se in(O) > 0, 


and we define the comparison functions by means of the 
equations 


p(t) =F Ye (0) +744(0) |, 


fax(t)=[ $ ve2(0) —- rp (0) | + 
+ 2re (0) | $ ra (0) — 5 Pa (0) |¢+-r4,(0), 


Finally, I’ will be considered to be the property that the 
inequalities 
Ir ,O— 2 OL< Fle, +7, 0) | 
ik ih al ih ik ’ 


hold, where Vin (0) = |r’, (0)]. 

To prove the inductiveness of I", we require the 
following lemma (G.A. Merman, 1953a). 

LEMMA 3.1. Jf the vector x(t) satisfies the inequality 


| r’—r | (0)<e 
in the interval OS t<7, then wt satisfies the inequality 


>> [v(0)—e]?2 +2r(0)[r/(0)—€] ¢+72(0), 
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for these same values of t, where v(0) =|r’(0)|. 
Proof: From the integral identity 


rr (0) +r’ (0) fie —r'(0)) dt 
we obtain 
F110) +r (OIL ft” —¥ ONaEl, 
Hence, by the hypothesis of the lemma we have 


r>|r(0)+r'(O)t]|—et, 


UF <a oP i 
or 


r2> [re (O) + r’ (0) t]? — 2e |r (0) + (O) ¢| e+ e2?, 
which proves the lemma when the inequality 
|r (0) +r’ (0) | <7 (0) + v (ODE. 
is applied. 
Using this lemma, we now show that if the property 
I’, i.e., the inequalities 
eA) <FPalO+re(O], (3.36) 


hold in the interval 0< <7, then the inequalities 


rie >| eta — Fre) | AF 


+ 2ra(O)[ Sr’, (0) — re (0) | é+72,(0). 8.37) 


will also hold for these same values of ¢. 
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To prove this fact, we apply Lemma 3.1 for which it 
suffices to set 


v= Vee, 
eee (04 (0) 4’, (ODI 


and to carry out some simple manipulations. 

From the expression (3.34) which defines the operators 
A;;, it immediately follows that if it is assumed that 
ri» (t)>f;,(t) in the interval O<t<7, then the inequalities 


‘ um] . ee | t f 
le dt ~ t d 
Nye \ > t ‘“ t ar + ‘\ my = mt + 
| rik mars | Pr es 7 2 < Mog fey 
0 J=90 6 j=) ‘ 
jris a Sk ” 


a n—1 : 
+d) tht Sim a, (3.38) 
0 


will hold for these same vaiues of ?. 

The results obtained show that the conditions of the 
second Theorem on the Inductiveness of Inequalities are 
satisfied by the case under consideration. Therefore, the 
inequalities (3.36), as well as the inequalities (3.37) which 
follow from them, are inductive properties for systems of 
gravitating bodies which satisfy the inequalities 


t n—] t F 

Liat 

sist (3.39) 
+ Sims] ez leaW)-+ra OV] <0. 


a 
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for all values of #20 

Moreover, (3.36) and (3.37) are automatically sat- 
isfied at the initial time ‘=0; therefore, a system of grav- 
itating bodies satisfying condition (3.39) will satisfy (3.36) 
and (3.37) for all values of +20. Thus, the distances be- 
tween the bodies tend to infinity as >. 

It remains for us to investigate conditions (3.39). 

All the integrals which appear in (3.39) are of the 
same type, namely, for any values of ¢ and 7 we have 
f { 
dt =| dt (3.40) 
() 


fag J asst? bast + ey? 
" 


where 


2 1 2 
ay =| Fry (0) — Fru (0) |, (3.41) 
by = [Fri (O)— zy () ]2ry (0), (3.42) 
Cy = Wy (0). (3.43) 
It is well known that the properties and the evaluation 
of any integral of the form (3.40) essentially depend on 
the sign of the discriminant 
Aj; = 4a, j€;;—53;. 
Using the values of a; ,, b; ;, and ¢;; given by (3.41) (3.43), 


as well as conditions (3.35), we find after some simple 
computations that 
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A, = 72, (0) [v3 (0) —r2 (0)] >0. (3.44) 


For A;;>Oand a,;;>0, the quadratic form a; j++; jt+ 
C;; IS positive for all values of t. Therefore, 


t 


dt de 
| ae 4- bet + cy <J Gg yt? + be gt + cay” (3.45) 


Furthermore, when A; ;>0 and, as in our case, a; ; >0, we 
have 
© 


j dt a 
J Ay + bette 


rt 2 b 
— wen are: stan a ‘: 


Vasy VAgy VOgy 


Since 6;;<0, as follows from equation (3.42) and condi- 
tions (3.35), from the last equation we obtain the estimate 


de < 2 
‘ Asjl? pe dy yl + ey VA; (3.46) 


Finally, from (3.40), (3.44), (3.45), and (3.46), we find 
that the inequalites 


<5 ou —. (3.47) 
e a Tuy (() Vv, (O) — rs (0) 
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hold for any values of ¢and jand for all ¢>0. From thi$ 
last result, we immediately obtain the estimate 


f F n—1 q n~1 P 5 
et. oN m, [at 
me | fer 7 Da, At a 7 ET < 
S#ik Jak 


r (VM —anig) (VE — mid 


— See (3.48) 
Fer (0) V0 (0) =r (0 


where M=)"%=§m;. 
On the basis of (3.48), we can state that if the 
inequalities 
r (M — m;) (M —m,) __ 


- Pe (0) VY v2, (0) — re (0) 
—F [POAT (0)] <0, 


hold at the initial time ¢=0, then (3.39) and all other 
consequences of it will hold for all values of ¢>0. 

To conclude our investigation, we can formulate the 
following theorem. 

THEOREM 3.6. If the initial conditions 


r, (0) — 27), (0) >, 
pa (0) — r (0) > 0, 


3V3n (AM —m,)( As — oe 
rik (V) 


(eq (0) +r, (O)] V 2%, (0) —7, (0) > 
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are satisfied at t=0 by a system of gravitating bodtes, then 
the distances r;, between the bodtes tend to infinity as to. 

We shall now make some comments apropos of the 
theorem just proved. 

The criterion we have obtained for the complete 
scattering of a system of bodies as t->= requires each of 
the velocities v;,(0) to be sufficiently great in comparison 
with both the reciprocal of the distance r;,(0) and the 
radial velocity 7’;,(0). The first requirement is intrinsic 
to the problem, whereas the second is a consequence of 
the mathematical reasoning used in the derivation of 
the criterion. 

The theorem imposes the fewest conditions on the 
initial state of the system when 7’ ;,(0)=0. We note, in 
passing, that it is just in this very case that the criterion of 
Theorem 3.4 fails to answer the basic question. When 
ry’ ,(0)=0, the first and second conditions of Theorem 3.6 
are automatically satished and the third simply becomes 


v5 «(M— M— 
vt, (0) >® x ( =a my) 


If the relative and radial velocities are equal, namely, 
if v,,(0)=7';,(0), then Theorem 3.6 fails to produce an 
answer. However, in this case, Theorem 3.4 is appli- 
cable and the criterion established there is the most suit- 
able one. Thus Theorems 3.4 and 3.6 complement one 
another. 


Finally, we note that the criterion obtained above may 
be regarded as a weakened variant of the criterion for 
hyperbolic motion formulated by G.A. Merman (see the 
reference 1953a) for the particular case of the three-body 
problem. We would like to remark on this variant in the 
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following connection. The criterion of Merman is not 
very transparent, and its form does not permit us to tell 
whether the set of initial conditions satisfying this criterion 
is not empty. In our case, this is obviously so. More- 
over, it is not clear under what conditions Merman’s 
criterion is applicable; however, they cannot differ essen- 
tially from those under which Theorem 3.6 holds. 

9. Let us now consider a system of gravitating bodies 
Po, Pi, ..., P,—1. We describe the motion of the system 
using Jacobi coordinates. Set p=p,—1, where p,—; de- 
notes the distance of P,_ from the center of mass of the 
remaining bodies. Let p denote the vector joining the 
center of mass of the system Po, P;,..., P,—2 to the body 
Pi-t- Then 


— a +. MM -|} + Ke 55224 ’ (3.49) 


where i, j, and kK are unit vectors directed along the 
£1, M1, and ¢, 1 axes, respectively. Differentiating 
(3.49) twice with respect to time, we obtain 


9” = itt +j jn + key 
or, by the equations of motion (1.7), 


pee — fi iY aU cay jth =}: (3.50) 


Now, using the relations 
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n—2 
| Mery)  OFt, n—1 
— a ar eamainaas ST 
Kat —0 tml ee 
n—2 
Mimy) Ors, n—1 
—_w a ae ————a ry 
O%n—} i=0 Ps, a—1 on.) 
s—2 
MM) — OFi nm 
Xu—| 4=0 rt n—] 1 


we can write (3.50) in the form 


“ 4 < MYM y—1 ° or; n—] FFs wat n——1] s, n—1 | 
ee a ee ee 
_ OMe -1 “eI 


Integrating this equation with respect to time from :=0 
to :>0, we obtain 


p —p' (0) = 
$ n—2 
or or Or. 
=—— 1 pith i 7 n—1 ane f, x—1 St k : n—1 dt, 
Ham ed re 05, 0%, _, Va) 


Using the fact that 


i Or ¢. s—1 + j “s n—1 tks vf —1 


ft 
Gal | 


<1, 


we find from this that 


ly —P OI < — Sa di. (3.51) 
ee 
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If it is assumed that no a-approach of P,—1 to the 
other bodies occurs for ¢>0 (the definition of aq-approach 
was given in Section 4 of this chapter), then 


rin—-1> ap. 


Using this to strengthen inequality (3.49), we then obtain 


an 
ly Ol SEIS (3.52) 
0 


where 


The integral inequality (3.52) is of the same type as those 
in the system of integral inequalities (3.33) considered in 
the preceding section. Applying the second Theorem on 
the Inductiveness of Inequalities to (3.52) and repeating 
the argument used to prove Theorem 3.6, we arrive at the 
following theorem. 

THEOREM 3.7. If the inittal conditions 


v (0) re 2p" (0) > G, 
v7 (0) —¢" (0) >0, 
, SAC ame vo eM 
Ir (0) ++ (0)] Vr? (0) — 3" (0) > Sy 
are satesfied at t=9 by a system of gravitating bodies, where 


v(t)=Yr tui te, 
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and if noa-approach of Py,—1 to the other bodies of the system 
occurs for any value of t20, then p(t)->@ asi-o. There- 
fore, the distance of P,—1 from each of the other bodies will 
become infinite as to. 

10. We shall say that a system of gravitating particles 
is completely dissipative as io if all the distances be- 
tween the particles increase indefinitely as to. 

THEOREM 3.8. If at the initial moment of teme i=0, the 
conditions 


p'(0)>0, 


8H 


eas’ 59) 


p’(0)> max {=o 


are satisfied by a system of gravitating bodies Po, Py,..., Py—s; 
where H is the total energy of the system, and if there occurs 
no a-approach of P,—1 to the other bodtes for any value of 
120, then the distance of P,,: from each of the other bodies 
increases indefinitely as to, but the system of bodies Po, 
P,, ..., Py—1 on the whole zs not completely dissipative. 

That the distance of P,—, from each of the other bodies 
in the system increases indefinitely as t+» follows im- 
mediately from Theorem 3.5a. 

We must now show that the system of bodies in quest- 
tion cannot be completely dissipative as i>. 

From the energy integral 


1 
1 ™~\! ' , , 
=3 (FAW +e) 


| 


we have 
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U 2% a1 (Ey + Nyt Sa 7 oe —4, 


and using the obvious inequality 


*n— as oa jet th 2 1a 2), 


we find that 


U Dz ret” (t)—H 


From Theorem 3.5a it follows that p’(#)>3 p’(0) for all 
values of 20. Therefore, 


Uu> = a (0)—H 


Hence by (3.53) we have 
U>0 


and, therefore, U does not approach zero. This means 
the system of gravitating bodies Po, Pi, ..., P,1 cannot 
be completely dissipative as i>. 

11. In proving some of the theorems, namely, Theorems 
3.5, 3.5a, 3.7, and 3.8, we assumed there occurred no 
a-approach of P,1 (whose motion we were studying) to 
the other bodies Po, P:, ..., P,—2. It seems to us that 
the presence of such a condition in the problem is due to 
its intrinsic nature. 

In fact, consider the aggregate of the gravitating bodies 
Po, Pi,..., P,—1, and suppose that P,—, is located inside 
the aggregate. If the initial velocity of P,_, is sufficiently 
great, and if while moving inside the aggregate P,,— does 
not approach the other bodies closely, then P,—1 will leave 
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the aggregate and go off to infinity. In other words, its 
distance from each of the other bodies will become in- 
finite as ->o&, 

However, if P,,; comes close to the other bodies, this 
could result in such an exchange of energy between the 
bodies in the system, that P,:, by having its velocity 
changed, could not leave the aggregate in spite of the 
fact its initial velocity was sufficient to have this happen 
under the other conditions. 

However, it is necessary to observe still one other 
aspect of the problem. It could happen that two of the 
bodies P; and P,—1 are ejected from the aggregate, and in 
moving an infinite distance away from the other bodies 
of the system, are within an a-proximity of one another 
or else arrive at this state after a certain time has elapsed. 
In this case, the theorems we have proved fail to solve 
our problem. 

Thus, the criteria formulated in Theorems 3.5, 3.5a, 
3.7, and 3.8 yield a solution to the problem only if the 
body in question experiences no subsequent a-approaches 
to the other bodies. Since these theorems do not indicate 
under what conditions any a-approach is excluded, they 
merely have qualitative significance. Nevertheless, they 
can be used quantitatively in some problems in the dy- 
namics of stellar systems. For example, by knowing the 
parameters of a stellar system, it is possible to estimate the 
probability of an a-approach of a given body to any of 
the other bodies and, therefore, to estimate the prob- 
ability of occurrence of the type of motions we have been 
considering, under given conditions on the initial state 
of the system. 

Finally, we note that the condition that there be no 
a@-approaches can be omitted in certain instances since it 
is automatically satisfied. Namely, this is so in the 
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application of Theorem 3.8 to the particular case of the 
three-body problem. It was in this way that the author 
first obtained criteria for hyperbolic-elliptic motions in 
the three-body problem. Later on, G. A. Merman, using 
our method of continuous induction, gave several inter- 
esting theorems containing more exact criteria for hyper- 
bolic-elliptic motions in the three-body problem. 


CHAPTER 4 


Method of Invariant Measure 


1. The application of the method of invariant measure 
to mechanical problems and, in particular, to the many- 
body problem is based on the phase interpretation of 
mechanical systems. 


The phase interpretation of a system of gravitating 
bodies is as follows. In Chapter 1, we saw that the 
second-order differential equations of motion can be re- 
placed by an appropriate number of first-order equations, 
the general form of which may be written as 


dz 
qo Xi (a, nis Be); 


ee ee ee ee ee (4.1) 


dz, 
ie »(%, a | 2) 


The state of the system of bodies is then completely 
defined at each moment of time ¢ by the values of vy real 
continuous functions of t, x1, x2,...,%, We-call each 
possible set of values of xi, x2, ..., x, a phase of the 
dynamical system. While regarding the collection of 
numbers x1, %2, ..., x, as ithe phase of the dynamical 
system, we shall, at the same time, think of it as a point in 
v-dimensional euclidean space £’, which we refer to as 
phase space. 

We now choose some moment of time as our initial 
time which we designate to be ¢=0. Each point pC’, 
with the exception of those which are singular points of 
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the given system, is a possible phase of the dynamical 
system at timet=0. The phase of the system at any time 
¢ other than the initial time will in general be different 
from the initial phase. We let f(p,t) denote the point of 
E that represents |this new phase, thus f(,0)=. As the 
time varies positively, f(p,t) will be defined in some in- 
stances for all positive values of #, whereas in ather 
instances, there will be a t>0O such that f(p,t) will be 
defined only for all positive t<r. A similar situation 
will prevail when the time varies negatively. 

As the time varies, the phase x1, %:,..., x, will change 
in accordance with the laws inherent in the motion of the 
mechanical system. Corresponding to these changes of 
phase, the point f(,¢) will move in a given way in E” and 
will describe a certain trajectory. 

By simultaneously considering all the points of phase 
space and their motions, we obtain the picture of the flow 
of the “phase-point fluid.”’ This picture describes the 
aggregate of possible motions of the dynamical system. 
By applying the method of invariant measure, we shall 
be able to characterize the motion of the particles in 
physical space in terms of the motion of the points of 
multidimensional phase space. 


2. For each pair of points of the space E” 


P= (%y Te, 00) My 
PY = (21, %y oy %) 


we introduce the nonnegative number 


(0, PI=N Ba? 
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called the distance between p’ and p’’. The general 
properties of the distance in £” are exactly the same as 
those in three-dimensional space: 

(1) (Identity property) p(p’,p’")=0 if and only if 
p’ _ p": 

(2) (Symmetry property) p(p’,p") = p(6",p’); 

(3) (Triangle property) For any three points p’,p’’ ,p 
of FE”, we have 


p(?’, p')+p( 0", b'”’\2> p(t’, bp’), 


Let » be a point and A an arbitrary set of Z’. The 
nonnegative number 


u(p, A)= infp(p, 9) 
qEA 


at 


is called the distance between the point and the set A. 

Using the notion of distance, it is easy to introduce into 
E” the notions of the neighborhood of a point and the 
neighborhood of a set. If p is any point of £’ and é€ is 
an arbitrary positive number, then the set of all points ~’ 
for which p(?,p’)<e is called a spherical neighborhood 
of # of radius € and is denoted by S(p,e). If A is any set 
of £’, then the set 


S (A, &) aay S (p, ®) 


is called a spherical neighborhood of A of radius e. 

Although it is impossible to give a graphical representa- 
tion of the motion of a point and its trajectory, or the 
motion of a set, in phase space, nevertheless, both the 
above interpretation with its geometrical terminology and 
the analogy of the properties of the space EF” to those of 
ordinary three-dimensional space considerably facilitate 
the discussion. 
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Since distances and neighborhoods are characterized in 
y-dimensional space by properties analogous to those of 
distances and neighborhoods in ordinary spaces, the theory 
of point sets and the theory of measure can be suitably 
generalized to the space EZ”. A presentation of these 
theories and, in particular, the basic concepts that we 
shall use may be found in books on the theory of sets and 
the theory of functions of a real variable. 


3. The general properties of the phase motions /(,t) 
will, of course, depend on the character of the functions 
X ;(%1,...,%,) that appear on the right-hand side of (4.1). 
For a wide class of equations of the form (4.1), to which 
the equations of motion of a system of gravitating bodies 
also belong, it can be shown that the phase motions have 
the following properties: 

(1) Uniqueness property—every point p of phase space 
(with the exception of singular points) defines a unique 
motion of f(,t) in an interval of time of one of the follow- 
ing four types: 


—o<dic+on, 

—o <i”, 
Uc<t<+to, 
Pct<t”. 


Thus, through each point of phase space passes only 
one phase trajectory. 

(2) Stattonarity property—in the interval in which 
f(p,t) is defined, the condition 


S(f(d, bh), #)=flb, +t); 


is satisfied. 
(3) Continuous dependence of the motion on time—for 
every ¢ in the interval of time in which f(p,t) is defined and 
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for any €>0, we can find a 6 >0 such that 


P(S(d, t), (Pt) <e 
for all ¢’ such that 
| t—t’ | <6; 


(4) Continuous dependence of the motion on the initial 
point—if the motion of f(p,2) is defined over the interval 
[41,42], then for € >O we can choose a positive 6 sufficiently 
small so that for any point g for which f(q,¢) is defined and 
satisfies the condition 


(f(a, 4), f(P, 4)) <6, 


the motion of f(,¢) exists in the interval [t,4] and satisfies 
there the condition 


pP(F(g, t), f(p, t))<e. 


4. Of great importance below will be the concept of the 
measure of a set of points of phasespace. We shall use the 
measure normally defined in multidimensional euclidean 
space by the direct generalization of the concept of 
measure in three-dimensional space. This measure is 
only meaningful for bounded sets. However, along with 
bounded sets, we are going to consider unbounded meas- 
surable sets both of finite and of infinite measure. To 
our end, it suffices to define the measure of an unbounded 
set as follows. Let Ri, Ro,...,R,,... bean unbounded 
monotonic increasing sequence of positive numbers. Let 
S, denote the sphere of radius R, with center at the origin. 
We call an unbounded set A measurable if each of the sets 
A()S,,ismeasurable. Thesequence meas(4()S,), »=1,2, 

. , either tends to a limit or becomes unbounded as 
n—>o. If the limit exists, we then define 
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meas A = lim meas (A()S,), 
m— & 


and in the contrary case, we set 
mes A = © 


Clearly, when this definition is applied to a bounded set, 
we obtain its usual measure. 

The great importance of the concept of measure is in 
the fact that Liouville’s Theorem on Invariant M easure 
is applicabie to the motion of a system of gravitating 
bodies in phase space. The invariance of measure im phase 
space is the basis that allows us to ascertain some of the 
very important properties of dynamical systems. 

THEOREM 4.1. (Liouville’s Theorem) If the equatzons 


az 
Fe ki (am, seey My); 


B68 Ss eee Yn es B (4.1) 


are such that 
aX, _ 
> == 


and if A is a measurable set of E” of finite measure sttch that 
for any point p © A the motion of f(p,t) ts defined for all 
t20(t<0), then 


meas f (A, t)=meas A 


(4.2) 


for all t>0(#<0), i.e., the measure 1s invariant when 
#20(#<0). 
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Proof: Let A be a measurable set of E’ of finite measure, 


where for any p€A, the motion of f(%,4) is defined for all 
t>0¢<0): 


meas { (A, t)= | dz, ... dx,. (4.3) 


f(A, 8) 
Consider the solution of (4.1) 
Ly 2 (t; a1, «22, 4), 
ae es de Sal a (4.4) 

Uy By (bi Ar, vey Or) 

satisfying the initial condition 
L —- a, 
for (=U, 


nl ¢ 


where (a1, ..., @,) is any point of the set A. 
Let us introduce new variables of integration m, ..., x, 
defined by equations analogous to (4.4): 
Vy = 2 (2; Zi, ea eg Ly), 


e ee @® oe (4.5) 


in which ¢ is regarded as a parameter. In terms of the new 
variables, (4.3) becomes 


meas {(A, t)= |D-dz, ... ah, (4.6) 


A 


where D is the Jacobian 
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DO v2 ++) Fa 
0 (2,, ao | 2.) : 


Differentiating (4.6) with respect to the parameter t, we 
obtain 


d | dD 
meas f(A, t) = \ar aa . . aby, 
A 


Then going back to the old variables x,,..., x, we find that 


+ meas} (A, = | Joe dn vee Qty. (&-T) 


(A ‘ f\ 


By the rule for differentiating a determinant, 


av > (4.3) 
Tt 7 2 


Where 


Ox, rt. 
ala vay Bievege TA ° 
? r) 


D _ Ty, ties oe Z.) 
q 


oF by 1) 


t,) e 
ten Xi te ‘.) 
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~~ OX 
Di Si Thy 
: > OL, ( 
¥ then noting that 


te for hes; 

é ok Zi+ik? on £y) bD 

eee ES 0 for kumi’ 
) =! 


Fy 4, Peay oe +) 


We fing 


Substituting this expression for D; in (4.8), we obtain 


av pV 
Ot) amy Ot, ' 


izpl 


and hence 


Using this last result, we can express formula (4.7) in the 
following form: 


(zal 


ev : 
A meas f(A, t) _ Dita + -te (69) 
f (4, #) 


This is the remarkable formula of Liouville. 
We now proceed to prove Liouville’s Theorem. From 
formula (4.9) and condition (4.2), we have 
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d = 
az meas {(A, t)=0, 


Integrating this with respect ¢ from ¢=0 up to ?>0 and 
taking into account that f(A,0)=A, we find that 


meas f (A, !)=meas A. 


The theorem is proved. 

It is necessary to emphasize the importance of the 
requirement that for any point €.A ,the motion of f(,t) is 
defined for all ¢>0(¢<0). If this requirement is not 
satisfied, the set f(A,t) and its measure may not be 
meaningful for all values of 220(¢<0). This can even be 
true for differential equations of a very simple form as the 
following examples show. 

EXAMPLE 1. Consider the motion defined by the 
equations 


dz __os re 4.10 
ae ; Wt ye ° ( ) 


Phase space for this system is the euclidean plane £?. 
From (4.10), we find that 


— —=—Yy, Y¥=Ce’, 


where C is the constant of integration. Thus the trajec- 
tories of the phase motions are exponential curves. Let 
us choose an arbitrary point p in the plane with coordinates 
Xo,yo, and let us consider the motion of f(/,t). Integrating 
the first of equations (4.10) subject to the initial condition: 
%=%9 when t=0, we find that 


x = —In(e—*o—t). 
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From this result it is easy to see that x40 as t->e—*o. 
Therefore, the phase point f(p,¢) goes off to infinity in a 
finite interval of time, and for all ¢2¢~*o, the motion of 
f(p,t) is not defined. This is true for any choice of the 
initial point p€ E?. 

EXAMPLE 2. Consider the motion defined by the 
equations 


(4.11) 


Phase space for this system is again the euclidean plane E. 
From (4.11), we find that 


dy 

ae 
ens 

yo=Ce ?, 


where C is the constant of integration. Thus the trajec- 
tories of the phase motions are normal curves. Let us 
choose an arbitrary point / in the plane with coordinates 
X0,Yo, and let us consider the motion of f(,t). Integrating 
the first of equations (4.11) subject to the initial condition: 
x=) when t=0, we find that 


t= \ 


From this we see that 
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to at 
X+»++0o as (—>—= [« 2 daz<cl, 
& 


Thus, the phase point f(p,¢) goes off to infinity in a finite 
interval of time when the time varies both positively and 
negatively. In this example, every motion is defined in 
an interval of the form 


<t<e’’, 
where ¢’ > —1 and ?’’ <1. 


5. The main object of study in general dynamics are 
the invariant sets of phase space. A set A is called postiwe 
(negative) invariant or invariant as t->o (t+->—o) if for 
every point p€A, f(p,t) is defined for all #20(¢<0), and 
for these times 


S(Pt)EA. 


A set A is called simply invariant if it is invariant both 
as t->o and ast->—o. It is clear that the union or in- 
tersection of any collection of sets that are invariant either 
as i> or as t->— ©, as well as the difference of two in- 
variant sets, is also invariant as © or as t>— ©. 

Let us consider the set of all points p€E for which 
the motion of f(p,t) is defined for all 420 (¢<0), and let 
= * (#7) denote this set. Obviously, 2+ (2-) is a 
maximal set, invariant as to (t>— ©); we call it the 
regular kernel of the dynamical system for >= (t->— ~ ), 

We call the set 
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cd ah (4.12) 


simply the regular kernel of the dynamical system. 
THEOREM 4.2. &t, #7, and = are sets of type Ga, 
1.e., sets which are the intersection of a countable number of 
open sets. 
We first prove the theorem for the set+. Let 


Ty Toye eg. Tag oe 


be an unbounded monotonic increasing sequence of pos- 
itive numbers as no. Let 7 denote the set of all 
points p of phase space for which the motion of f(p,¢) is 
defined in the interval [0,7,]. If p€&, then due to the 
continuous dependence of the phase motion on the initial 
point, we may choose a positive number 6, such that for 
every point g, with g€S(p,5,), the motion of f(q,t) is 
defined for all z¢€[0,7,]. However, this means that 
S(q, 6 »\ Ee ,and therefore = is an open set. Obviously, 


on 


belt] = + oe 
Ce ie fe 


from which it follows that “+ is a set of type G; and is 
therefore measurable. 
The theorem is proved in an analogous way for =. 
The validity of the theorem for the set = follows 
immediately from (4.12). 


6. Let us take some set A of EH” at the initial moment 
of time. Each point p€A will move in its own way in E’, 
and at some other moment of time, all the points p of A 
will occupy new positions f(~,/) and form a new set; we 
denote this set by f(4,?). 

Let p and A be such that f(#,¢) and f(4,t) are defined 
for all#20(¢<0). From the definition of f(p,t) and f(A,t). 
it follows’ that if p€A, then f(p,t)€f(4,t) for every 
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t>0(¢<0). Conversely, if f(p,¢)Ef(A,t) for some 
t20(t<0), then pEA. 

From this remark, we immediately have the following 
simple facts. If Ai and Ag are sets such that A:()A2.=9, 
and if f(A1,¢) and f(As,t) are defined for all :>0(¢< 0), 
then for these values of ¢ we have 


FAN f(A2t) = 0. 


On the other hand, if 4; and Azare sets such that 41(\42= 
0, and if f(Axt) and f(A2,t) are defined for all :>0(¢< 0) 
then for these values of ¢ we have 


f(A, OM (Ase, #)40. 


A set A is called recurrent as to (t->— © ) if for any 
T>0, we can find a ?’>T such that 


Af\f(A,t')4#0 (ANf(A,—t’)#9). 


Conversely, the set A is not recurrent as ->o> ([>—©) 
if we can find a 7 >0 such that 


ANf(A, H=0 (ANS(A,—4)=0) 


for all i>r. 

Theorem 4.3. If a set A ts recurrent as [> (t->— © ), 
and tf f(p,t) ts defined for all t<O(t>0) and for any potnt 
DEA, then A ts recurrent as t>— © (t+). 

If A ts not recurrent as t--o(—>— @), and tf f(p,t) ts 
defined for all t<O(t>0) and for any point pCA, then 
A ts not recurrent as t->— © (ito). 

Proof: Let the set A be recurrent as +>. Then, for 
any positive 7 no matter how large, we can finda !>T 


such that 
AQ f(A, ') #0. 
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We then have 


L(A,—t) OSHA, #),—#) #0, 
f(A, —U)NA +0, 


and this implies that A is recurrent as t>— , 

For a set A which is recurrent as t->— ~, the theorem 
is proved analogously. 

The second part of the theorem is easily shown by 
arguing by contradiction. 

THEOREM 4.4. Ifa set A ts such that f(A,t) is defined 
for all real values of t, and uf we can find a number +t >0 
such that A(\f(A,t)=0 or AC\f(A,—t)=0 for all t>7, 
then the sets 


Arun f(A,—nr), 7 ~»J(A,—7T),A JS(A,T), os f(A, nT), eects 
are patrutse disjoint. 
We prove the theorem by contradiction, and we assume 
that we can find two integers 7 and j, with 7<j, such that 


F(A, ir)Q f(A, jr) 40. (4.13) 


Consider the sets A and f(A,(j-2)T); since 7-121, we have 
by assumption 


At) f(A,(i—t)r) =0. 
and, therefore, 
flA, ir) Nf (f(A (JT), 17) =, 
However, 


F(f(A ’ (j—t)T), iT) =f(A, jt). 
Thus 


fA unin f(A, j7) = 
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and this contradicts (4.13), and so the theorem is proved. 
The theorem is proved analogously for the case where 
Af) f(A,—t)=0 for all t27. 
7. A point p is called recurrent as to (t->— © ) if for 
any positive T no matter how large and any e>0, there 
exists a time t’ > T(t/<—T) such that 


f(b, EVE S(p. €). 


A point p is called simply recurrent if it is recurrent both 
when to and when t>—-. 

THEOREM 4.5. If the point p is recurrent as t->© 
(t»—), then for any h>O(t<0), the point f(p,ti) ts 
also recurrent as t--> © (t->— © ), 

We shall only prove the theorem for a point p which is 
recurrent as i>. The case where the point is recurrent 
as t>—-o is treated analogously. 

Let us take arbitrary positive numbers € and 7, and 
let us choose 5 = 6(€,7) sufficiently small so that 


AS(p. 6) AVCSHA(P, 4),€). 


This is possible because of the continuous dependence of 
the motion on the initial point. Since the point p is 
recurrent as t->, we can find a ?’>T such that 


I (pt) €5(2, 8). 
If gis a point of the set A, then by the definition of f(g,t) 
and f(A,t) it follows that f(g,t)€f(A,t). Using this fact, 
we Can write 


SF. #), b) E F(S(P, 6). ACSF (D4), €), 


and since 


Sf (p. t’), A) =f (f(b. 4). 8), 


it follows that 
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fC (2, t), LES CG (p, &), €). 


This means that the point f(p,f1:) is recurrent as to. 
We thus sce that if # is recurrent as t->~, then all points 
f(p,t) will be recurrent as ¢-+o. Therefore it is natural 
to call the motion of f(%,¢) recurrent as t->o. The same 
remarks are also valid for points that are recurrent as 
t-+— oo. 

8. Consider the set P of all recurrent points in the phase 
space of a mechanical system. For the following, it will 
be important to establish its measurability. We shall 
also need to consider the set P+ (P~) of all points of phase 
space that are recurrent as t->© (or +— ©). 


THEOREM 4.6. P+, P-, and P are sets of type G; and 
hence are measurable. 


Proof: Let 


Fy Vay os 4: Tages s 


bean unbounded monotonic increasing sequence of positive 
numbers and 


Une Noa, --+) Nay - 


a monotonic <lecreasing sequence of positive numbers 
which tends to zero as a limit. Let Pt denote the set of 
points ~ for which one can find a T>7T, such that 
p(p, f(p.7))<n,. Consider a point p’€ Pt, and let 
p(p’. f(p’,7))=B <7,. We then choose a positive number 
€ <(n,—8)/2 and a 5>0 sufficiently small so that for 
any point q for which p(p’,q)<6, the inequality 
p(f(p’ .t),f(9,t))<€ is satisfied for all positive <7. This 
is possible because of the continuous dependence of the 
phase motion on the initial point. We then have 
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(9, D') KBE < "Jo (% — 9) 
o(p’, f(P,*)) =B, 


o(f (D’; t)y ACE <8 < aA (%,— B), 
and hence by the Triangle Property, it follows that 


0(9, 4 (91 ™)) < Tne 


This means that from the condition p’€ P+, it follows that 


S(p’,s) C Py, 


i.e., Pf is an open set. Obviously 
= + 
P= 1) Pi, 
comp] 


and this implies that Pt is a set of type Gs. 

The fact that P~ is also a set of type G; is proved 
analogously. 

The validity of the theorem for the set P then follows 
from the obvious relation 


P=Ptf\P-. 
Coro.uary. The sets E7—-P and =—P are measurable. 


9. The theorems presented above were of a preliminary 
nature. We now proceed to prove the main theorems. 

THEOREM 4.7. Let I be a measurable set not containing 
potnts that are recurrent as t-»—@ ), where for each PEM, 
f(p,t) ts defined for all t20 (or 1<0). If meas M>0O, then 
one can find a bounded set m*Cot of postive measure whtch 
4s not recurrent as t-> © (o7 t-»>—@ ). 

We shall only prove the theorem for ~~. The case 
t->—@ is treated analogously. 
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Let 


T1, Ta, ee» Tye yee 


be an unbounded monotonic increasing sequence of pos- 
itive numbers, and 


M1, Nar + + + May oe 


a monotonic decreasing sequence of positive numbers 
tending to the limit zero. We let SW, denote the points 
PEON such that 


F(D, EVE SUD, Mn) 


for all t27,. Since the set 9% does not contain any 
recurrent points, each point pC belongs to at least one 
of the 91, and, therefore, 


M = UM,. (4.14) 


Let us show that 91, is measurable. To do this, we shall 
prove that each of the sets Si, is closed with respect to IM. 
Let p* be a point of IN which is a limit point for M,. 
Let us assume that p*€OM,; then we can find a time 
i >T, such that 


p(p*, f(p*, t)) =8<y,. 


Now let €< (7,—8)/2, and let us chooseé >0 sufficiently 
small so that for any point g for which p(p*,q)< 6, the 
inequality 


p(f(p*, ti), f(g", i)<e, 
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is satisfied. This is possible because of the continuous 
dependence of the phase motion on the initial point, 
Moreover, the point g may be chosen in IN, since p* is a 
limit point of 9%,,. Therefore, 


0G P*)<b<e<F (m— 8), 
o(p*, f(p*, 4,)) =8, 


o(f(p*, 4), F(a A)<e <> (Nn — 8): 
or by the Triangle Property 


P(g, (4G, &)) <u, 


This contradicts the fact that g is in ,. Therefore, 
p*EOM, and I, is closed. 
From (4.14), it follows that 


> meas I, >> meas PR > 0, 


nut 
and, therefore, we can find a set IN, such that 
meas Qi, >. 


Let x be a point of metric density of the set I,, and 
consider the set 


Mi, =WeNS (x, Fre) - 


It is clear that Sj is measurable and that meas 91%, >0. 
The diameter of 91% does not exceed 9,/2. But for our 
point p of 314, we have 
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pp, f (BP, t)>) np 
for all t>7,. Therefore, 


MAAR = 


for {>7T,. Setting 9t*=9N; completes the proof of the 
theorem. 

THEOREM 4.8. If A ts a measurable set in EY such that 
any subset of tt of positive measure 1s recurrent as t->o 
(or t->— © ), then all points of A, with the possible exception 
of a set of measure zero, are recurrent as t->@ (or i~>—@ ). 

We prove the theorem by contradiction, and we assume 
that there exists a set SW of positive measure all of whose 
points are not recurrent as t>© (or t>—o). Then by 
Theorem 4.7, we can find a bounded measurable set 91* 
of positive measure which is not recurrent as t>o 
(or +>—o ). However, this contradicts the conditions 
of the theorem. 

Let A be a measurable but otherwise arbitrary set of 
phase space of the dynamical system. We call the set A 
regular as t->o if for any 7 >0, we can find a set A* such 
that meas A* < ~ and 


U f(A,t)C A. 
9gtg 


A set which is regular as —» is defined in a similar 
way. 

Let B bea set and # a point of phase space; we call the 
point p recurrent with respect to the set B as i> (t>— ~) 
if for any T>0, we can find a positive (or negative) 
number ?#’, with (¢/|>T7, such that f(p,t’)€B. 

THEOREM 4.9. Let A be a regular set of points of phase 
space with invariant measure and W a set whose points are 
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recurrent with respect to A as t->o (or t-»—om). Zhen 
any measurable non-recurrent set M*E W ts a set of measure 
zero. 

We prove the theorem by contradiction, and we assume 
that there exists a set I7*C W which is not recurrent and 
which has positive measure. Then we can give a positive 
number 7 such that the sets 


f(M*, nt), n=0,1,... 


are pairwise disjoint. Since A is regular, there exists a 
set A*, with meas A*< ©, such that 


U f(A, t)C A®. 
<tr 


0 


It is not difficult to see that for any pEIN*, we can 
choose a subsequence from the sequence f(p,nT), n= O,1,2, 

. , all of whose terms are points in A’*. 

Let us consider the sets 


Me = f(M", nt), 


a=0,1,2... 
and let 
D, = A° 1M. 
Each D, is measurable, being the intersection of two 


measurable sets. Moreover, the sets D, are disjoint. 
Therefore the series 


>) measD, (4.15) 
v=0 


converges since 
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D, Cc A*, 


fee] 
U DC A*, 


m0 


meas U D, <measA* << o 
and - 
Symeas’ =measl) D,< a, (4.16) 
Further, 
f (Dy, — nt) = M*N /(A*, — az) (4.17) 
is the set of points p in IN* for which we have 


f(p, nt) A®. 
Each point of 31* belongs to an infinite number of the sets 


(4.17) and, therefore, to each of the sets 
f+ °) 
Rs =U (D,, igs | (4.18) 
nwi1,%, ... 


Thus, for any n we have 


M* Cc H,. (4.19) 
By the invariance of measure, we obtain 
meas {(D|, —vt)=meas D (4.20) 


and so by (4.18) and (4.20) 


meas 9%, == Dj meas/(D,, —vt)= >) meas D. 


Now, the series (4.15) converges and 
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is its remainder. Thus for any €>0, we can choose an 
integer N= N(e) sufficiently large that 


meas in < ? (4.21) 


for all 22 N. 
Taking into account relation (4.19), and inequality 
(4.21), we conclude that 


meas IN* < :, 


and since € is an arbitrary positive number, this implies 
that 


meas W* =). 


We have arrived at a contradiction, and this proves the 
theorem. 


10. We now study the motion of a system of m grav- 
itating bodies in phase space making use of the analytical 
relations of celestial mechanics and the above general 
measure theoretical results on dynamical systems. 

Let us first consider the set Q of all points P€ & for 
which 


max {7;;} > © as [> —O, 
ij 


max {7,;}-+» « as f—+»to. 
‘fj 
Weare not going to investigate the measure of Q. For our 


purposes, it will suffice to prove that Q satisfies the 
following purely logical theorem. 
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THEOREM 4.10. From the assumption that 


meas Q>0, 


at follows that all the points of Q, with the possible exception 
of a set of measure zero, satisfy the condttion 


inf (min {7,;}) = 0. 
t>0 sf 
Proof: Let Q- denote the set of all points p¢ E™ for 
which 
r®*(t)= max {r,j}— 00 as t>—o, (4.22) 


and let us show that this set is measurable. From condi- 
tion (4.22) it follows that for any point p€Q7 and any 
number R>0O, we can find a positive number T so large 
that 


r¥ (t)2R for all t <—T. 
Let 

Ray Pay oie le os 

Die Digs eee, Dig S88 


be two unbounded monotonic increasing sequences of 
positive numbers. We let Q~{Rz, Tm} denote the set of 
all points for which 


r* (t) DR, for all t<—T, 


We now show that this set is closed. 
Consider the set G, of all points #€ E™ for which 


r*(p) Oe 


since r*(p) is a continuous function of the phase point, 
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G, is an open set. Now arguing by contradiction, we 
assume that O~{Rz,Tm} is not closed. Then we can find 
a point p*€Q-{Rz,7,,} which is a limit point of this set. 
For this point, we can find a moment of time ¢’/<—7 
such that 


f(p*, tyeG, 


Since the set G, is open, it Is possible to find a number 
€>0 sufficiently small that 


S(f(p*, t'), #) CG. 


Due to the continuous dependence of the phase motion on 
the initial point, we can choose a 6 >0 sufficiently small 
that 


f(S (p*, 2), LY CS (f(p*, C), €) Gy. (4.23) 


However, ~* is a limit point of the set] Q-{R,,7,,} and 
therefore we can find a point p’ such that 


p é Q- (Ry, Las 
p'ES (p*, 4). 

From (4.23) and (4.24), it follows that 
f(D t’) E Gy, 


which leads to a contradiction. Thus the set Q7{R,,T7,,} 
is closed. 
Let us next consider the set 


(4.24) 


Qr= U (@- (Ry Tn) 2). (4.25) 


Its measurability follows from the results obtained above. 
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Now, Qj is the set of points of phase space each of which 
eventually leaves the set 


i < ft, (4.26) 


permanently as t-->—o., 
From the obvious formula 


= 1 65 


it follows that the set Q~ is measurable. The set 07 con- 
tains the points of phase space that each eventually leaves 
any set of the form 
r®* (P) <u (4.27) 
Te 3 ee Tee 

permanently as i>—o. 

Carrying out a similar argument for t»«, we can 
construct measurable sets 


OF= U (Cr {Res Tn) 2) 
analogous to the sets (4.25). Then 


Q- N(=- Or) 


is a measurable set consisting of the points in phase space 
each of which eventually leaves any set of the form (4.27) 
permanently as +>+—« and which is also recurrent with 
respect to the set (4.26) as ->o. 

It is further clear that 


Q= U (2-NE-@)) 
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from which it follows first of all that Q is measurable. 
We are assuming meas Q>0, and therefore, we can find 
a v such that for the set 


Q,=Q-/N (H—-QF) (4.28) 
we have 
meas Q,> 0. 
It is very important to note for the following that Q, is 


the set of points each of which, as t>—o, eventually 
leaves any set of the form 


r* (pb) © Ras 


wel, 2, 


permanently, and as ¢+->@ is recurrent with respect to 
the set 
r® (P) SR. 


Let 


Qi, A2,--+, Any --e 


be a monotonic sequence of positive numbers which con- 
verges to zero, and let A, denote the set of points p in 
E®, for which 


ty (p) =min {r;;} 2a, for all #20. 


Let us show that each A, is closed. Consider the set Gj of 
all points for which 


Te (b)< ar: 
since ry(p) is a continuous function of the phase point, 


Gis open. Let us assume that A, is not closed. Then 
we can find a point py € A, which is also a limit point of 
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thisset. For this point, we can then find a time ¢’’>0 such 
that 


f (Darr #°) © G'e. 


Since Gj is an open set, we can find an €>0 sufficiently 
small such that 


Sf (Par YY) CGp- 


By the continuous dependence of the phase motion on the 
initial point we can choose a 6>0 so small that 


F(S (Par 8), OD OS (CF (Dar #1), €) CG. (4.29) 


Now, x is a limit point of A,. Therefore, we can find 
a point p” such that 


p”’ E Ag, 


4.30 
Pp’ ES (dx, 5), ( 


From (4.29) and (4.30) it follows that 
fe", VEG, 


and this leads to a contradiction. The set A, is closed 
and hence measurable. 
We now consider the proof of the theorem itself. 
We argue by contradiction. If the theorem is not 
true, then we can find a p» such that for the measurable set 


=, A,; 


we have 


meas @ > 0. (4.31) 
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Let us take an unbounded monotonic sequence of positive 
numbers 

i, Ha, eee neo cues 
and let us denote by W, the set of all points whose phase 
coordinates satisfy the inequality 


{ a ; 4 | 
|#\=\> S m (ar + yi, +2 = Ut ity. 
i==l 
Since the left-hand side of this inequality is a function of 
the phase and continuous on #, each of the W, is an open 
set. Writing Q) as the sum of a denumerable number of 
measurable gets 


a= 0 (en,) 


and taking into account inequality (4.31), we conclude 
that we can find a positive integer h such that for the set 


W= Qe a) W,, 
we have 
meas W >0. 
For the following, it will be important to enumerate 
the properties of the points ~ in the set W: 
(1) if pE W, then as t+—», f(p,t) leaves any set of 
the form 
r* (pr) SR, 


nazJ, 2... 
permanently, i.e., r*->(t)~ as b>—o ; 
(2) every point p€ W is recurrent with respect to the set 


r* (p)<R, 
as [>o; 
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(3) if p€W, then the motion of f(p,t) is such that 
ry (t)2 Oy for all #20; 

(4) if p€W, then the motion of f(p,t) satisfies the 
inequality 


ww 
1 ; 2 19 
|| =|a >) m (22+ y?4 2) Ul <M, 


>| 


for all real values of ¢. 
We let B denote the set of points p€ & for which 
min {744} 2 Ay for all 20. (4.32) 
|| < H,; (4.33) 


From the preceding considerations, it follows that this set 
is measurable. Moreover, it is clear that if p€ W, then 


f(p, t) € B for all #20 (4.34) 


so that B is invariant as t>«. Let us show that the 
phase velocity v, considered as a function of the points of 
phase space, is bounded on the set B. The square of this 
velocity is given by 


oS HES) a+) 4) 


7 = ies sei) tras ays) +(e see) tar tate, | penne 


From inequality (4.32), it follows that the potential func. 
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tion U is bounded from above; therefore we can find an 
H*>0 such that U< H*, and hence we have 


" 
4 c , 12 
= Sin, (22+ yo +z. )\<H+H"*, 
=) 
i.e., the kinetic energy is bounded from above. Thus the 
functions x2, 2, and 2? are bounded from above. Dif- 
ferentiating U with respect to x;, we obtain 


oU > mymy = Oris 
Se aa OO 2 S ames : ’ 


Since 


aul <i 


we have that 


ise |S yee 


From this inequality and condition (4.32), it follows that 
av \2 
(a) 


is bounded from above. In a similar way, we can show that 
the functions 

(— 2 (= )" 

wun) > (dae) 


are bounded. From these results and equution (4.35), it 
follows that v?, and therefore also v, is bounded on the set B&B. 


Let 
A= BS (nk), 
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where S(”R,) is an open sphere with center at the origin, 
being a positive integer (we are considering the motion of 
the 2 bodies in a coordinate system with origin at the 
center of mass). We have shown that the phase velocity 
is bounded on the set B. Therefore, we can find a positive 
number b such that |v|<b. Since A€B, this estimate is 
also valid for the set A. Let us choose a T>0 but other- 
wise arbitrary; for any p€B, the length of the path tra- 
versed by the point} f(p,!)/in phase space in the interval of 
time from 0 to 7 does not exceed 67. Thus every point 
PCA is such that for all ¢ in the interval 0< ¢< T, we have 


f (p,t) € S(mR, +57). 
It then follows that 
A, t)C S(nR,-++ dr), 
Uf (As 2). S(AR, +b) 


and since S(nR, +57) is a bounded set, meas S(nR, +07) 
< oo and therefore, A is regular as [->~™. 

To complete the proof of the theorem, we must obtain 
a contradiction. 

Let us first show that there exists a set 90*C W which 
is not recurrent as !>—-~ such that meas IN* > 0. Assume 
that this is not true. Then every measurable subset of 
W of positive measure is recurrent. However, in this case, 
we have by Theorem 4.8 that all points of W, with the 
possible exception of a subset of measure zero, are recurrent 
as i>—~. For the motion of f(p,¢), determined by any 
point DEW, we have that 


max {r,;}-»> co as t—-»—oo. 


This is impossible for recurrent points, and we arrive at a 
contradiction. Thus, there exists a set M*CW, with 
meas 90*>0, which is not recurrent as t»—o. By 
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Theorem 4.3, SIU* is also not recurrent as t->o. Now, 
we know that each point p€ W is recurrent as [>© with 
respect to the set 


r*(p)<R,. 


Therefore, if at some moment of time the point f(,t’) 
belongs to the set 7*(p)<.R,, the distances between the 
gravitating bodies will not exceed R,. Therefore, the 
distance of any of the bodies from the origin, located at 
the center of mass, does not exceed 7, , i.e., f(p,t’)€ S(mR,). 
Hence, every point p that is recurrent with respect to the 
set r*(p)<R,, will be recurrent with respect to S(”R,). 
These remarks, in conjunction with condition (4.34), allow 
us to say that every point p€ W will be recurrent as i> 
with respect to the set 


A= BS (nf), 


a set that is regular as to. We have thus arrived at 
the conditions of Theorem 4.9 from which it follows that 
meas 9%*=0.|Our contradiction is obtained, and the 
theorem is proved. 

We shall call a system of gravitating bodies almost 
stable as to if 


ss {753} +0 as t->+o, 
‘ 


and if we can find a number a >0 such that the condition 
min {r.j} >a 
4 


is satisfied for all t2 0. 

We shall call a system of ” gravitating bodies stable 
as t->o if wecan find positive numbers Rand a such that 
for all +20 

ssh {ra} <R, = (ry) pe 
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At present, it is not known whether there exist almost 
stable systems which at the same time are not stable. 

We shall say that a system of m gravitating bodies 
contains a body which has come from infinity if 


é 
max {rij} 0co as t-—-+—o. 
i 


THEOREM 4.11. For all inttial condtitons, with the pos- 
sible exception of a set of initial conditions of measure zero, 
a system of gravitating bodies which contains a body that has 
come from infinity cannot be almost stable. 

We prove the theorem by contradiction. 

If the theorem is not true, we can find a number a >0 
such that the set Q* of all points p€@ satisfying the 
conditions 


max {ryjy}a~o as t—->—o, 
ij 
max {ryy}> oOo as ¢—-»- 0 
iy 


inf (min {r,,}) 22 
ij 


#>=0 


(4.36) 
is a set of positive measure: 
meas Q* > 0. (4.37) 


However, Q*CQ, where Q is the set of points pC x 
satisfying conditions (4.35). From the fact that Q*€Q 
and the condition (4.37), it follows that 


meas Q >0O. (4.38) 


Now, by Theorem 4.10, we know that when condition 
(4.38) is satisfied by all points of Q, with the possible 
exception of a subset of measure zero, the condition 
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inf (max {r,,)) = 0. 
(30 = éf 


must be satisfied. 
This result, clearly, is also valid for the subset Qe 
d this contradicts condition (4.36). > 


an The theorem is thus proved. 


CHAPTER 5 


Analysis of Some Cases of the Evolution 
of a System of Gravitating Bodies 


1. The two-century-old development of the planetary 
cosmogony has often been described as a haphazard 
alternation of diverse hypotheses, as a process of trial and 
error that has produced no definite results. This, how- 
ever, is not so. The development of cosmogony as a 
science is a process that has followed the general pattern 
of growth of knowledge. Cosmogonical knowledge has 
been acquired in the form of relative truths, and as these 
truths have evolved, we observe the gradual formation 
of the elements of an objective science, which thus be- 
comes continually more perfect. Of course, the growth 
of cosmogony, like that of other sciences, has been a very 
complicated and contradictory thing, in which true and 
false principles have emerged side by side. However, 
in the course of its development, certain ideas have 
appeared that have persevered, and with modification and 
expansion, have remained as part of the science in one form 
or another. The history of science shows that such evolv- 
ing ideas constitute that part of the science which, to 
the extent of its growth, form the elements of objective 
knowledge. 

In the cosmogony of the planets, one such idea is, 
first of all, that the solar system and naturally formed 
systems similar to it evolved out of rotating matter 
scattered in space. Much controversy arose about the 
physical state of the matter, about its genetical relation- 
ship to the various other cosmical bodies, and about the 
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nature of the process by which the matter was accumulated 
and transformed into large bodies. However, a general 
basis has stood the test of time and is being confirmed 
by the experience of the two-century-old development of 
the planetary cosmogony. This basis thus comprises the 
elements of objective knowledge. The idea that the 
planets were formed from rotating matter scattered in 
space follows almost inevitably from the properties of 
the solar system and the natural laws governing it. 

Thus the basic questions we have to consider in cos- 
mogonical investigations are the process by which the 
matter scattered throughout wide regions of space was 
accumulated in a small volume, the process by which this 
matter was condensed, and the process by which tiny 
particles were unified into large bodies when the cosmical 
systems were formed. There is no doubt that gravita- 
tional forces play an essential role in these phenomena. 
However, today, it is regarded as well established that 
gravitational forces alone are insufficient to explain their 
occurrence, and of fundamental importance also are phys- 
ical phenomena involving energy conversion. Therefore, 
the full analysis of these cosmogonical processes is beyond 
the scope of gravitational mechanics. 

Nevertheless, the theorems of celestial mechanics can 
be of considerable help in these investigations. We begin 
by considering the following. The rotating matter scat- 
tered in space, from which the solar system was formed 
and which we take as its initial state, i.e., as the object 
whose evolution is impending, can be regarded as a system 
of many gravitating bodies. The final product of evolu- 
tion, the solar system, is also a system of many gravitating 
bodies which are fewer in number and larger in size. 
Without going deeply into the processes and mechanism 
by which the initial system was transformed into the final 
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one, but by comparing them only as two gravitating 
systems of bodies, we shall discover that this transforma- 
tion cannot be realized solely on the basis of the laws of 
gravitational mechanics. By looking only at the mechan- 
ical side of the question, we cannot, to be sure, explain the 
physical factors which are important to the transformation 
under consideration. Nevertheless, a consideration of the 
mechanical side of the question does yield significant 
results. We may not be able to discover the nature of 
these physical factors, but we can get some idea about 
them by determining the mechanical consequences of these 
factors. 


2. Let us consider gravitating bodies Po, Pi,... , 
P,—1, and let r;; denote the distance between P; and P,. 
We shall say that a completely hyperbolic approach of the 
bodies occurs if all the distances satisfy the condition 


ri3(t)>+ © as t>—o 


It is not difficult to see that a completely hyperbolic 
approach is a type of motion that is actually possible for a 
system of ” gravitating bodies. Namely, we know that 
in purely dynamical problems, it is legitimate to reverse 
the sign of the time. Therefore, if the inequalities (3.24) 
hold for the case of time varying negatively, then from the 
appropriate reformulation of Theorem 3.4, it will follow 
that a completely hyperbolic approach takes place. More- 
over, since 7;; and r;; are continuous functions of the 
coordinates of phase space of the system of gravitating 
bodies, the set of points of phase space whose coordinates 
satisfy condition (3.24) is an open set. Hence, it follows 
that the set of points of phase space which are initial 
points for the “regime’’ of a completely hyperbolic ap- 
proach is a set of positive measure. 
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Let us now show that the motion in which the conditions 


1;;(t) > © as t>—o (5.1) 
and 
r,Q<R as t>0, (5.2) 


are satisfied simultaneously, where R is any suitable 
positive number, is impossible for a system of n gravitating 
bodies. 


Disregarding a set of initial states of measure zero (the 
case H=0), we can assert that from (5.1) and Theorem 1.1, 
it follows that H>0. Therefore, by Theorem 1.2 at least 
one of the distances between the bodies tends to infinity 
as t->«, which contradicts condition (5.2). 


We thus see that for a completely hyperbolic approach 
of gravitating bodies, a system cannot arise in which the 
mutual distances remain bounded for all >0. However, 
if the bodies approach each other, then the unification of 
not all, but some portion of them into a subsystem with 
bounded 7;; is still possible. This was rigorously proved 
for the three-body problem in the theory of capture de- 
veloped by O. Yu. Schmidt, G. F. Khilmi, N. N. Parisky, 
and G, A. Merman. For a larger number of bodies, it 
should be possible to give examples in which such motions 
occur, however, they would be somewhat artificial and 
would entail very formidable computations. At any rate, 
it is possible to get a sufficiently clear picture of our 
phenomena by means of the following considerations. 

If from a number of bodies that have come from 
infinity, a stable subsystem is isolated following the ap- 
proach of the bodies, then on the one hand, U+O ag 
(~>»—o, and on the other hand, there are distances be- 
tween the bodies that are bounded on thesett>0. There- 
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fore, we can find a constant U) and a time T such that 
U> Up for t>T. 


Then from the formula 
t ty 


P(t) = IR + 1% t+2Mf [U+ 2H) dit, 
vr OU 


in which we must take H>0, we see that following the 
approach of the bodies, the moment of inertia J?(¢) will 
increase more rapidly as > than as t-»—© despite the 
fact that all the distances increase indefinitely as i>—o, 
and only a portion of them increase indefinitely as to. 
This means that the gravitational association of a portion 
of the bodies is accompanied by the rise of motions 
involving increased velocities. 

The physical interpretation of these facts is as follows. 
The too high kinetic energy which the bodies have as the 
result of approaching from infinity prevents their complete 
gravitational association. However, when certain condi- 
tions are satisfied, a portion of the bodies can form a.stable 
subsystem. This can happen when the gravitational in- 
teraction of the bodies as they approach each other results 
in a redistribution of the energy causing some parts of 
the system to begin to move rapidly with respect to the 
others. If the kinetic energy of the relative motion of 
parts of the system constitutes a large percentage of the 
total energy, these parts can then become stable systems. 
The unification of a portion of the bodies into one or 
several stable subsystems is necessarily accompanied by 
the ‘‘draining off’’ of the energy in the system, i.e., by 
the appearance of relative motions of the subsystems in 
which the surplus of energy preventing any association of 
the bodies is spent. Theorem 3.8 shows that if at least one 
of the bodies in a system of ” gravitating bodies has 
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sufficiently high kinetic energy, the other bodies cannot 
comprise a dissipative system. 

We thus see that the association and scattering of 
gravitating systems are connected in the closest way. 
If one of these phenomena occurs, it will be compensated 
for by the occurrence of the other. Association and 
scattering in gravitating systems are thus divergent aspects 
of the one phenomenon of the gravitational interaction of 
matter. 

The above results do not mean that whenever we have 
a hyperbolic approach of many bodies, it is impossible 
for all of them to be unified into a stable system under 
their mutual attraction. Our investigations merely point 
out that this is impossible as a purely dynamical process. 
However, if we make allowances for processes of a more 
general type, then this may conceivably happen. 

In fact, if whenever the bodies drawing near under 
their mutual attraction gives rise to or is accompanied by 
some physical process in which the surplus of kinetic 
energy preventing their association is absorbed, then 
association can be realized without the phenomenon of 
scattering or else with a relaxation of it. In this case, 
scattering is ‘‘replaced”’ by a physical process that absorbs 
the surplus of mechanical energy and converts it into 
non-mechanical form. 

For meteoric dust-clouds, the most important form of 
“drainage” of mechanical energy occurs when a portion 
of it is converted into thermal energy during non-elastic 
collisions of the particles. When there issucha ‘‘drainage,”’ 
association is conceivable without the scattering of dust 
particles or with the relaxation of scattering. In par- 
ticular, it is possible for meteoric dust particles to be 
captured when a star passes a meteoric dust-cloud closely. 
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The dynamics of stellar passages through a cloud of 
meteoric dust particles was investigated by T. A. Agekyan 
who worked out a theory of capture tied up to the particles’ 
loss of kinetic energy in mutual nonelastic collisions. 


3. We can formulate some further general results by 
using Theorem 4.11 which was proved in Chapter 4. To 
do this, we must know under what circumstances the 
evolution of any particular system of real bodies is 
describable by means of the dynamical laws for a system of 
n gravitating bodies. 


Real bodies such as the stars, planets, small solid 
particles, etc., are idealized as mass points. However, the 
idealization of real bodies as mass points is permissible 
as long as the distances between them are sufficiently great 
in comparison to their dimensions. Once the bodies are 
distributed so closely that the distances between them are 
comparable to their dimensions this idealization ceases 
to be correct. In this case, various physical phenomena 
arise which are accompanied by the conversion of mechan- 
ical energy into a non-mechanical form. Finally, when 
approaching closely, these bodies may undergo considerable 
changes: they may collide or be disintegrated or, con- 
versely, they may be unified, etc. Therefore, whenever 
the motion of real bodies is to be described in a given 
interval of time by the dynamical equations for a system 
of gravitating bodies, we must assume that the distances 
between the bodies do not become smaller than some 
positive number a whose value will depend on the concrete 
properties of the given system of bodies. 


With these remarks, let us consider a system of grav- 
itating bodies Po, Pi, ..., P,—1, which we assume to have 
the following properties: 
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oj(t) > co as L->—oO, 


=1.2,....0—] 
rey(t) +> Cc as t—»-+-oo, 
4, J 1,2, 2.0) wm 
$31 


i.e., we consider a system of bodies P;, Ps, ..., Po 
which does not get scattered as +> © , to which is adjoined 
the body Py coming in from infinity. We next show that 
if a physical system is formed after Po is adjoined to the 
system consisting of Pi, Ps, ..., P,~1 which is stable as 
t->o, then the process by which this system is formed 
could not be purely mechanical: it is necessarily accom- 
panied by the conversion of mechanical energy into non- 
mechanical form. 

In fact, Theorem 4.11 implies that except for a set of 
initial conditions of measure zero, for every system of 
many bodies containing a body incoming from infinity 
but none outgoing to infinity, the greatest lower bound of 
the mutual distances r;; is zero. Therefore, arbitrarily 
close approaches of the bodies in the system will occur, 
which means also the nonmechanical phenomena that 
accompany them. 

For the sake of simplicity, we have considered the case 
where the system contains only one body coming in from 
infinity. However, our conclusions still hold for any 
number of such bodies. 

4. We now derive an inequality which we shall use 
below. 

We consider # gravitating bodies P;, Ps,...,P, with 
masses 1, M2, ...,m,, and we describe their motion 
with respect to a set of coordinates with origin at the 
center of mass of the system. 

Consider the following identities whose validity can 
easily be verified directly: 


EVOLUTION OF A SYSTEM OF GRAVITATING BODIES 105 


aby ba 
= 5 [ea — 92+ Wa — tad + a, — 22) +r 


éax 1,2, ot 


Multiplying these identities Sosa by m1, mo,... 
and m,, and summing, we obtain 


27=P+0?, (5.3) 


where 7 is the kinetic energy and P and (2? are given by 


V2 Padi G7 Y,— Y, 2) + 


=a 


cas a (VA — tH, > y4 <j (2 — 82) (5.4) 


«=! 


= mz. (5.5) 
é=) 


Now, consider the first sum in the expression on the 
right-hand side of (5.4). The variables x;, y;, x, and y; 
appearing in thissum are not independent. They are con- 
nected by relations (1.6) expressing the area integrals. If 
we use Lagrange multipliers to determine the relative mini- 
mum of the first sum in equation (5.4), we find that it is 
equal to c?/J?. 

In an analogous way, we find that the relative minima 
of the two other sums on the right-hand side of (5.4) are 
respectively 


cand 


yo Ot 
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Using these results and recalling that 
Pad+a+a4, 
we find that 
Q? 


P2 ih (5.6) 


From relations (5.3) and (5.6) it follows that 


2T 2 yl +?, 
or 
rey 
JI?> 5.7 
> (5.7) 
From the energy integral, we have that 
2T = 2(U+H); 


substituting this expression for 27 in (5.7), we obtain the 
required inequality 


pp —_% (5.8) 
2(U+ A) —Q? 

5. A system of a large number of gravitating bodies is 
called stationary, if despite the movement of its component 
bodies, the distribution function for the distances between 
the bodies remains practically constant and in this sense 
is an invariant of the motion. Clearly, a system can be 
stationary only if its total energy is negative. 

If a system is stationary, then 

2 J 2 


J+=const; at =(Q, 
at? 
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Therefore, it follows from the Lagrange-Jacobi equation 
that 


U+2H=0. 
From this equation and the energy integral, we have that 
T=—HA, U=—2H, 


or taking into account that H <0, we can write 
T=|A|, U=2|H|. (5.9) 


6. Let a physical system S be in the state (A) at time 
to in which it is a collection of m) mutually gravitating 
bodies. Suppose that at time 41>, the system S passes 
into a new state (B) in which it is a system of m, gravitating 
bodies. As to the transition of S from the state (4) to 
the state (B), we assume nothing other than it has occurred 
as the result of the gravitational interaction of the moving 
bodies. Of course, of essential importance in this transi- 
tion are also the nonmechanical phenomena which are not 
manifested in the dynamical laws for gravitating bodies. 
Under this formulation of the problem, the number of 
bodies m) and m, in the respective states (A) and (B) are 
not necessarily equal. 


We let 
To. Uo, \ol, Jé, 88, Q5 
and 
T,, U,, |, J?, 87, 2? 
denote the values of the quantities 7,U, |JI|, J?, Q?, and 


©? for the respective states (A) and (B) (in other words, 
at the times f and ¢;). 


Suppose that state (A) satishes the condition 


108 QUALITATIVE METHODS IN THE MANY BODY PROBILEM 


avjt 
Tr as! (5.10) 


and suppose that after transition into state (B), S is a 
stationary system of bodies occupying a region so small 
that 


&% 


ee ee (5.11) 
2 |Ho| — QF 


Fix 


Let us study the transition of S from state (A) to state 
(B). From the Lagrange-Jacobi equation (1.14) and con- 
dition (5.10), it follows that 


ly, +- 2/1, <0, 
From this inequality and the energy integral, we obtain 
—H,, a / igre 


Taking into account that JJ <0, we have that 


\H,| = 7, 
or that 


2) Ho) ZT. 


Subtracting 02 from the left and right-hand sides of this 
inequality, we obtain 


2 (Hol — 25 > 7,,—2i 
Now, from the definition of 7) and (22, it follows that 


and therefore that 


EVOLUTION OF A SYSTEM OF GRAVITATING BODIES 109 


2|H.|— 23 >0. 
Thus if 6%>0, the following inequality is satisfied : 


@ 
5.12 
iH a>” sie 


This inequality shows that condition (5.10) is math- 
ematically consistent (admissible). 

With the system S in state (B), we now write inequality 
(5.8) using the condition 


(7, = 2|H,!. 
which holds hecause of the stationarity of state (B). We 
obtain 
ey 


? 1 
We Dh | at (5.13) 


From inequalities (5.11) and (5.13), it follows that 
2|H 1-2? ~ 2]My| — 22 


Now, in order for this inequality to hold, at least one of 
the following inequalities must be satisfied: 


|\H,|>lAgl, (5.14) 
Bi< 8}. (5.15) 
Qi < Ob. (5.16) 


The results obtained, of course, do not permit us to 
clarify the character of the non-mechanical processes that 
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occur in the transition of system S from state (A) to state 
(B). However, some information concerning these proc- 
esses can, nevertheless, be obtained. Namely, it is pos- 
sible to indicate some of the mechanical consequences of 
the physical phenomena that accompany the transition 
of S from state (A) to state (B). 

Thus, let us consider inequalities (5.14), (5.15), and 
(5.16). If inequality (5.14) is satisfied when S passes 
from state (A) into state (B), this means that the transi- 
tion is accompanied by phenomena that cause a decrease 
in the mechanical energy of the system. Thus, for ex- 
ample, if the system in the state (A) is a cloud of solid 
particles, then the nonelastic collisions between these 
particles will cause a decrease in the mechanical energy 
of the system since as the result of these collisions a portion 
of the kinetic energy is converted into heat. 

If, however, inequality (5.15) is satisfied, this means 
that the transition from state (A) into state (B) is accom- 
panied by phenomena for which the moment of orbital 
momentum of the bodies comprising the system S in the 
state (A) passes into the moment of axial rotation of the 
bodies which make up SS in the state (B). The transition 
of the orbital moment into the axial moment is possible 
if there are collisions between the bodies, or if many 
bodies are unified into a single larger body. 

We must still consider the case where inequality (5.16) 
is satisfied. If all the orbits of the bodies in the system 
were circular, (2=0. As the orbits deviate more from 
circles, i.e., as the radial velocities r’ become larger, the 
value of 922? will become larger. 

Therefore, we can regard (2? as a measure of the 
deviation of the actual orbits from circular orbits. But 
then the fact that (5.16) holds, means that the transition 
of S from state (A) to state (B) is accompanied by a 
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rounding off of the orbits, i.e., by the orbits changing into 
a more circular shape. This can occur as the result of the 
interactions and disintegrations of the bodies when they 
collide, or when many bodies are unified into a single 
larger body. 

It is understood, of course, that one cannot exclude the 
possibility that the phenomena may be so related that 
after the transition of the system S from state (A) to state 
(B), all three inequalities (5.14), (5.15), and (5.16) will be 
satisfied. 

The analysis of the problem in question by means of 
the methods of celestial mechanics does not permit us to 
obtain more concrete results. To this end, it is necessary 
to carry out further cosmogonical investigations. How- 
ever, we must point out that good agreement does exist 
between our results and those given by present-day 
cosmogony. The investigations of academician 0. Yu. 
Schmidt, confirmed later on by the work of L. E. Gurevich, 
A. I. Lebedinskii and other research workers, have shown 
that the formation of the planetary system from scattered 
matter rotating about the sun was accompanied by con- 
versions of the kinetic energy of this matter into thermal 
energy, as well as by the unification of the scattered matter 
into the planets, and by the rounding out of their orbits. 
The conversion of mechanical energy into nonmechanical 
form turned out to be the basic process, and the two other 
processes its consequences. 
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